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ABSTRACT 


Data  foe  twenty-seven  geomorphic  and  coastal-process 
attributes  were  collected  at  l-km  intervals  for  800  kilometers 
of  the  mid-Atlantic  barrier  coast  between  Cape  Benlopen, 
Delaware,  and  the  North  Carolina-South  Carolina  border. 
Correlation  and  principal  component  analysis  was  run  on  fifteen 
of  these  attributes  in  an  attempt  to  classify  the  coast. 

Local  subregions  (between  55  km  and  190  km  in  length) 
showed  organization  and  interrelationships.  These 
relationships  are  not  as  clear  when  the  entire  800-km  data  set 
is  considered  in  the  same  analysis,  indicating  that  coastal 
geomorphic  and  process  systems  are  in  adjustment  to  local 
environmental  conditions  to  a  greater  extent  than  they  are  to 
regional  conditions. 

The  large  number  of  variables  resulted  in  a  classification 
of  the  mid-Atlantic  coast  into  twenty-four  distinct  barrier 
types  based  on  process  and  morphology.  A  coarser 
classification  of  the  area  identifies  seven  types  based  on 
attributes  of  coastal  strike,  sediment  size,  offshore  slope, 
wave  frequency,  shoreline  erosion,  inlet  frequency,  and 
offshore  bars. 


LIST  OF  TABLES 


LIST  OF  FIGURES 


Figure  1.  Map  of  the  eastern  United  States  showing  the  location 
of  the  800  km  study  area  between  Cape  Henlopen, 
Delaware,  and  the  North  Carolina-South  Carolina 
border. 

Figure  2.  Index  map  of  the  study  area  showing  the  six 
geographic-geomorphic  subregions  discussed  in  the 
text.  Also  shown  are  the  thirty  segments  used  as 
orientation  segments  used  in  analysis  of  wave,  wind, 
and  shelf -si ope  data. 

Figure  3.  Strike-parallel  plots  of  shoreline  strike,  island 
width,  lagoon  width,  overwash  penetration  distance, 
and  rate  of  shoreline  eroson.  For  details  of  data, 
see  Appendix  A. 

Figure  4.  Strike-parallel  plots  of  tidal  range,  storm  surge, 
inlet  frequency,  frequency  of  3-m  dunes,  and  sediment 
size.  For  details  of  data  see.  Appendix  A. 

Figure  5.  Strike-parallel  plots  of  the  frequency  of  1-m  waves, 
the  frequency  of  3-m  waves,  offshore  slope  to  5-m 
depth,  offshore  slope  to  9-m • depth,  and  mean  bar 
number.  For  details  of  data,  see  Appendix  A. 

Figure  6.  Morphology  of  longshore  bars  along  the  mid-Atlantic 
coast. 

Figure  7.  Schematic  classification  of  the  mid-Atlantic  barrier 
coast  based  on  principal  components  analysis  of  the  15 
variables  in  Table  4.  Each  illustration  represents 
the  mean  state  for  that  segment  of  the  coast. 
Although  the  sketches  do  not  conform  to  actual  scale, 
they  do  show  relative  scale  between  the  segments. 
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INTRODUCTION 

In  1971  we  began  an  investigation  of  regional-scale 
variations  of  the  sedimentary  landforms  along  the  Atlantic 
coast  under  sponsorship  of  the  Office  of  Naval  Research.  To 
date  we  have  reported  on  variations  in  coastal  landforms  (Dolan 
et  al.  1975);  offshore  bathymetry  (Resio  et  al.  1977);  barrier 
island  topography  (Vincent  et  al.  1976) ;  inshore  bathymetry 
(Dolan  et  al.  1977);  equilibrium  profiles  (Felder  et  al.  1979) , 
coastal  marine  fauna  (Hayden  and  Dolan  1976) ,  and  shoreline 
erosion  and  shoreline  configuration  (Dolan  et  al.  1977;  Hayden 
and  Dolan  1979;  Dolan  et  al.  1979) ,  and  more  recently,  Atlantic 
coast  wave  climates  (May  et  al.  1983). 

The  substantial  data  inventories  developed  from  these 
studies  provided  the  basis  for  classifications  of 
regional-scale  coastal  environments  and  landform  types  (Dolan 
et  al.  1975),  and  barrier  islands,  lagoons,  and  marshes, 
(Hayden  and  Dolan  1979) .  The  resulting  classification  units 
were  on  the  order  100  km  to  300  km  along  the  coast. 

More  recent  data  collections  on  coastal  processes  and 
responses  were  designed  to  analyze  regional-scale  associations 
at  higher  resolutions.  Accordingly,  an  800-kilometer  section 
of  the  Atlantic  coast  was  investigated.  Fifteen  attributes 
spaced  at  1-km  intervals  form  the  data  matrix  which  we  analyzed 
using  numerical  classification  procedures.  The  results  are 
reported  here. 
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MID-ATLANTIC  MICROTIDAL  BARRIER  COAST  CLASSIFICATION 


DATA  COLLECTION  METHODS 

Xhfi  Study  Aifi.a 

The  coastline  under  study  extends  approximately  800  km  from 
Little  River  Inlet,  North  Carolina/South  Carolina,  to  Cape 
Henlopen,  Delaware,  and  encompasses  the  mid-Atlantic  barrier 
islands  of  North  Carolina,  Virginia,  Maryland,  and  Delaware 
(Figure  1,  Table  1).  Within  the  study  area,  800  sample  sites 
were  designated  at  1.0-km  intervals.  The  sites  were  numbered  1 
to  800  moving  from  south  to  north,  and  were  identified  by 
coordinates  of  latitude  and  longitude.  Each  site  was  also 
assigned  a  map  and  transect  number  corresponding  to  the 
University  of  Virginia  Orthogonal  Grid  Mapping  System 
(O.G.M.S.) .  The  O.G.M.S.  map  and  transect  number  specify  the 
location  of  each  site  to  the  nearest  100  m  along  the  coast  on 
base  maps  prepared  from  U.S.  Geological  Survey  7-1/2  minute 
(1:24,000)  series  topographic  maps.  The  most  recent  maps 
available  were  updated  with  1976  aerial  photography  to  show  the 
shoreline  position  on  a  common  date. 

Selection  of  Variables 

The  first  phase  of  this  project  was  the  identification  of 
quantifiable  attributes  of  the  barrier  coastline  of  North 
Carolina,  Virginia,  Maryland,  and  Delaware.  The  principal 
constraints  on  selection  of  the  variables  were:  (1)  variables 
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Figure  1.  Map  o£  the  eastern 
United  States  showing  location 
of  the  800  k&  study  area  between 
Cape  Henlopen,  Delaware,  and  the 
North  Carolina-South  Carolina 
border.  . — ■ 
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TABLE  1 


MID-ATLANTIC  BARRIER  ISLANDS  AND  BEACHES 

Cape  Henlopen,  DE 
Rehoboth  Beach,  DE 
Fenwick  Island,  DE-MD 
Assateague  Island,  MD-VA 
Wallops  Island,  VA 
Assavoman  Island,  VA 
Metomkin  Island,  VA 
Cedar  Island,  VA 
Parramore  Island,  VA 
Hog  Island,  VA 
Cobb  Island,  VA 
Wreck  Island,  VA 
Ship  Shoal  Island,  VA 
Myrtle  Island,  VA 
Smith  Island,  VA 
Fishermans  Island,  VA 
Cape  Henry,  VA 
Virginia  Beach,  VA 
Sand  Bridge,  VA 
False  Cape,  VA 
Currituck  Banks,  NC 
Bodie  Island,  NC 
Pea  Island,  NC 
Hatteras  Island,  NC 
Ocracoke  Island,  NC 
Portsmouth  Island,  NC 
Core  Banks,  NC 
Shackleford  Banks,  NC 
Bogue  Banks,  NC 
Hammock  Island,  NC 
Browns  Island,  NC 
Onslow  Beach,  NC 
Ashe  Island,  NC 
No  Name  Island,  NC 
Figure  Eight  Island,  NC 
Shell  Island,  NC 
Masonboro  Island,  NC 
Carolina  Beach  Island,  NC 
Smith  Island,  NC 
Oak  Island,  NC 
Holden  Beach  Island,  NC 
Hales  Beach  Island,  NC 
8unset  Beach  Island,  NC 
Bird  Island,  NC 
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had  to  be  quantifiable  to  be  suitable  for  statistical  analysis; 
(2)  data  on  each  variable  for  all  segments  of  the  800-km  study 
area  had  to  be  accessible  through  maps,  bathymetric  charts, 
aerial  photographs,  existing  University  of  Virginia  coastal  data 
sets,  or  published  literature;  and  (3)  values  for  each  variable 
had  to  be  assignable  to  sites  at  1-km  intervals  throughout  the 
study  area. 

Twenty-seven  variables  were  identified  for  data  collection 
and  analysis  (Table  2)  .  Eight  of  the  variables  represent 
physical  processes  acting  in  the  coastal  zone,  and  the  other  19 
variables  are  geomorphological  attributes  of  the  barrier 
coastline  and  adjacent  water  bodies. 

Resolution  of-Varlables 

The  twenty-seven  variables  in  the  coastal  classification 
system  can  be  divided  into  two  distinct  subsets  based  on  the 
spatial  resolution  of  the  data  that  corresponds  to  each 
variable.  Eighteen  of  the  variables  are  classed  as  high 
resolution  geomorphological  variables.  The  other  nine 
variables,  eight  of  which  are  physical  process  attributes,  are 
considered  low  resolution  variables.  The  only  low  resolution 
geomorphic  attribute  is  the  slope  of  the  continental  shelf. 

The  arbitrary  distinction  between  high  and  low  resolution 
was  made  with  the  guideline  that  high  resolution  variables  have 
been  systematically  measured  at  1-km  intervals  along  the  coast, 
and  their  values  are  typically  not  constant  over  long  stretches 
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of  coastline.  Quantifiable  variation  is  usually  present  in  high 
resolution  variables  along  a  10  to  20  km  length  of  coastline. 
For  example,  barrier  island  width  can  be  measured  at  1-km 
intervals  on  1:24 ,000-scale  maps,  and  typically  varies  along  10 
km  of  shoreline.  In  contrast,  the  frequency  of  1.5-m  high  waves 
is  considered  a  low  resolution  variable  because  it  varies 
slightly  along  small  stretches  of  ocean  shoreline,  and  local 
wave  data  is  not  available  for  much  of  the  800  km  study  area. 

Data  Variables 

Geomorphic  Variables  (High  Resolution) 

gtogllne  string 

The  strike,  or  orientation,  of  the  shoreline  was  measured 
at  each  1-km  interval  site.  A  visual  "best-fit"  straight  line 
was  drawn  over  the  kilometer  of  shoreline  south  of  each  site  on 
1:24, 000-scale  U.S.G.S.  topographic  mapB.  The  strike  of  this 
line  was  measured  in  degrees  east  of  north.  As  examples,  a 
north-south  striking  coast  facing  east  (i.e.,  the  ocean  is  to 
the  east  of  the  barrier  island)  has  a  strike  of  0°;  a  coastline 
striking  east-west  and  facing  south  has  a  strike  of  90°;  a 
north-south  striking  coast  facing  west  has  a  strike  of  180°;  and 
an  east-west  striking  coast  facing  north  has  a  strike  of  270°. 
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Inlet  Frequency 


The  along-the-coast  frequency  of  inlets  was  measured  at 
each  site  by  counting  the  number  of  inlets  within  a  24-km 
segment  of  coastline  centered  on  the  site.  The  mean  length  of 
mid-Atlantic  barrier  islands  is  approximately  12  kmi  therefore 
24-km  segments  were  used  for  inlet  frequency  counts  to  reflect 
the  wide  variability  in  inlet  frequency  along  the  mid-Atlantic 
coastline.  The  inlets  were  counted  on  l:24r000-scale  D.S.G.S. 
topographic  maps  updated  with  air  photos  to  show  the  1976 
shoreline. 

TflBggiaghY 

Five  different  measures  of  barrier  topography  were  made  on 
the  most  recent  available  U.S.G.S.  l:24,000-scale  topographic 
maps:  (1)  frequency  of  dunes  between  1.5-2. 9  mr  (2)  frequency 
of  dunes  between  3. 0-4. 5  m,  (3)  frequency  of  dunes  between 
4. 6-6.0  m,  (4)  frequency  of  dunes  higher  than  6.0  m,  and  (5) 
frequency  of  dunes  higher  than  3.0  m.  These  values  were 
determined  at  each  site  by  examination  of  the  topography  of  a 
1.0-km  segment  of  the  barrier  to  the  south  of  the  site.  A 
transparent  grid  was  placed  over  the  l:24r000-scale  base  map, 
dividing  the  kilometer  under  study  into  ten  100-m  segments 
stretching  across  the  island  normal  to  the  shoreline.  The 
highest  topographic  contour  intersecting  each  of  the  ten 
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transects  was  recorded,  and  these  data  were  converted  into 
values  for  the  five  topographic  intervals  listed  above.  Only 
the  most  seaward  dune  field  was  considered  in  the  counts. 
Back-barrier  features  such  as  Jockey's  Ridge,  North  Carolina, 
were  not  considered. 

As  an  illustration  of  the  procedure,  if  the  maximum 
elevation  at  five  of  the  ten  transects  for  a  particular  site  was 
3.0  m,  in  three  transects  it  was  1.5  m,  and  in  two  transects  it 
was  less  than  1.5  m,  then  the  values  for  the  five  variables 
would  be  as  follows: 

(1)  1. 5-2-9  m  dune  frequency  «  30% 

(2)  3. 0-4. 5  m  dune  frequency  *  50% 

(3)  4. 6-6.0  m  dune  frequency  *  0% 

(4)  greater  than  6.0  m  dune  frequency  »  0% 

(5)  greater  than  or  equal  to  3.0  m 
dune  frequency  ■  50% 

Using  this  method,  the  range  of  possible  values  for  each 
variable  in  each  case  is  0%  to  100%. 

Offshore  Slope 

Offshore  slopes  to  5.5-m  and  9.1-m  depths  were  determined 
from  the  most  current  editions  (1981  or  1982)  of  National  Ocean 
Survey  1:80 ,000-scale  bathymetric  charts.  The  water  depth  was 
divided  by  the  horizontal  distance  from  the  shoreline  to  the 
appropriate  bathymetric  contour  measured  normal  to  the  shoreline 
strike  to  obtain  the  slope  in  m/km. 
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Bar  Number 


Thirteen  sets  of  color  infrared  aerial  imagery  of  the 
mid-Atlantic  coast,  between  1970-1979,  were  examined  to 
determine  predominant  longshore  bar  patterns.  None  of  the  photo 
sets  covered  the  entire  study  area.  On  the  average,  there  were 
four  photo  sets  used  for  bar  analysis  at  each  site.  The 
plan-view  bar  morphology,  as  evidenced  by  the  number  and 
position  of  breaker  lines,  was  mapped  from  all  photo  sets. 
Shoreline  attachment  points  of  the  bars  were  also  noted.  After 
photo  interpretation  was  completed,  the  mean  number  of  bars  and 
the  standard  deviation  of  the  bar  number  was  calculated  for  each 
site.  (The  complete  bar  analysis  will  be  described  i.n  a 
separate  ONR  technical  report) . 

Rate,  of  Shoreline-Change 

The  mean  rate  of  shoreline  erosion  or  accretion  at  each 
site  was  measured  by  the  Orthogonal  Grid  Mapping  System 
(O.G.M.S. ) .  This  system  involves  enlarging  air  photos  of  the 
coastline  to  a  common  scale  (1:5,000),  matching  landmarks  on  the 
enlarged  photos  with  features  on  1:5,000  base  maps  prepared  from 
U.S.G.S.  topographic  quadrangles,  and  tracing  the  shoreline 
(mean  high  water  line)  and  vegetation  line  onto  a  transparent 
overlay.  The  mean  and  standard  deviation  of  the  rate  of 
shoreline  change  is  calculated  by  measuring  the  distance  between 
an  arbitrary  base  line  on  the  base  maps  and  shorelines  traced 


from  photos  taken  on  various  dates.  Thorough  explanations  of 
the  O.G.M.S.  procedure  and  its  reliability  are  given  by  Dolan  et 
al.  (1978,  1980). 

For  sites  within  the  study  area,  the  rate  of  shoreline 
change  was  determined  from  up  to  seven  sets  of  aerial  photos, 
dated  193  8  to  1980.  At  all  sites,  photos  spanning  at  least  25 
years  were  mapped.  In  the  O.G.M.S.,  shoreline  positions  are 
digitized  at  100-m  intervals  along  the  shore,  and  rates  of 
change  are  calculated  at  each  of  these  100-m  transects.  To  make 
this  data  compatible  to  the  1-km  interval  sites  used  in  this 
project,  the  O.G.M.S.  shoreline  change  values  at  the  ten  100-m 
transects  immediately  south  of  each  site  were  averaged  to  yield 
a  mean  and  standard  deviation  of  the  shoreline  rate  of  change 
for  the  site. 

Qvfiryflgh.  £engtifl.tlgn .  Distance 

The  distance  between  the  shoreline  (mean  high  water  line) 
and  the  vegetation  line  was  also  measured  on  the  photos  used  in 
the  O.G.M.S.  This  distance  is  known  as  the  overwash  penetration 
distance  (OPDX) ,  and  it  was  calculated  for  each  of  the  1-km 
interval  sites.  As  with  the  rate  of  shoreline  change,  values 
measured  at  100-m  interval  transects  were  averaged  over  1-km 
stretches  adjacent  to  each  site  to  produce  a  value  for  the  mean 
and  standard  deviation  of  the  overwash  penetration  distance  at 
each  site. 
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Overwash  penetration  distance  has  also  been  called  the 
active  sand  zone  width  or  unvegetated  beach  width.  "Active" 
refers  to  the  dynamic  nature  of  the  unvegetated  beach.  On  this 
part  of  a  barrier,  aeolian  sediment  transport  is  a  continual 
process,  and  overwash  is  an  important,  though  infrequent, 
process.  Many  interrelated  coastal  attributes  must  be 
considered  in  site-level  or  regional-scale  interpretation  of 
overwash  penetration  distance  data,  including  the  frequency  and 
character  of  overwash,  the  rate  at  which  washover  deposits  are 
revegetated,  dune  stabilization,  engineering  structures,  and 
other  human  alteration  of  barrier  islands. 

Sediment  Size 

A  lengthy  review  of  the  literature  on  mid-Atlantic  beach 
sediments  revealed  that  there  was  no  single  high-resolution 
sediment  sampling  study  that  covered  the  entire  800-km  study 
area.  Consequently,  average  sediment  size  data  were  extracted 
from  four  sources  in  the  literature  and  from  the  unpublished 
results  of  four  coastal  field  studies.  Table  3  lists  the 
sources  of  sediment  data  used  for  each  segment  of  the  study 
area. 

The  sediment  data  sources  were  chosen  on  the  basis  of  the 
extent  of  the  sampling,  the  spacing  of  the  sample  sites,  the 
part  of  the  beach  that  was  sampled,  and  the  nature  of  the  grain 
size  statistics  that  were  calculated.  In  all  cases,  the  samples 
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were  collected  from  the  upper  foreshore  near  the  berm  crest. 

There  are  numerous  beach  sediment  studies  of  single  sites 
or  individual  islands  reported  in  the  literature  that  provide 
detailed  information  for  small  segments  of  the  coast.  However, 
for  the  purposes  of  this  study  we  chose  the  data  sources  that 
are  most  likely  to  reveal  the  regional-scale  variability  of 
beach  sediment  textural  trends  along  the  mid-Atlantic  coast, 
while  keeping  sampling  biases  to  a  minimum. 

Island Width 

Island  width  was  measured  at  each  site  on  1 : 24 ,000-scale 
U.S.G.S.  topographic  maps  that  had  been  updated  with  air  photos 
to  show  the  1976  shoreline  configuration.  The  width  from  the- 
ocean  shoreline  to  the  bay  shoreline  was  measured  perpendicular 
to  the  shoreline  strike. 

For  mainland-attached  beaches,  such  as  Virginia  Beach, 
Virginia,  the  mean  overwash  penetration  distance  was  used  as  a 
substitute  measure  for  island  width.  The  rationale  for  this 
substitution  is  that  the  mor phodynamic  character  of  the 
unvegetated  portion  of  a  mainland-attached  beach  is  functionally 
similar  in  many  ways  to  that  of  a  barrier  island. 

frftgppn  Width 

The  width  of  the  water  body  (lagoon,  bay,  sound,  etc.) 
separating  the  barrier  from  the  mainland  was  measured  at  each 
site  on  1  s 24 , 000-scale  U.S.G.S.  topographic  maps.  These 
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distances  were  measured  from  the  bay  shoreline  of  the  barrier  to 
the  bay  shoreline  on  the  mainland.  The  measurement  at  each  site 
was  made  perpendicular  to  the  ocean  shoreline  strike  of  a  3-km 
segment  centered  on  the  site.  For  mainland-attached  beaches  the 
lagoon  width  is  zero. 

Shelf  Slope 

The  slope  of  the  continental  shelf  was  determined  by 
measuring  the  distance  from  the  ocean  shoreline  to  the  183-m 
(100-fathom)  depth  contour.  The  measurements  were  made  on 
small-scale  National  Ocean  Survey  nautical  charts  (1:416,994  and 
1:1,200,000  scale).  Shelf  slope  is  a  relatively  low-resolution 
attribute  compared  to  the  other  geomorphic  variables  because  it 
does  not  vary  appreciably  between  adjacent  1-km  sites.  Rather 
than  measuring  shelf  slope  at  each  of  the  800  study  sites,  the 
study  area  was  divided  into  30  "orientation  segments”  by  drawing 
best-fit  lines  on  a  1:1, 000, 000-scale  map  of  the  region  (Figure 
2).  Within  each  orientation  segmei.c  the  coastline  is  nearly 
straight.  The  endpoints  of  these  segments  are  points  that 
coincide  with  major  shifts  in  shoreline  strike. 

The  distance  from  the  shoreline  to  the  183-*  depth  contour 
was  measured  along  a  line  perpendicular  to,  and  bisecting,  each 
of  the  30  orientation  segments.  The  vertical  difference  (183  m) 
was  then  divided  by  the  horizontal  distance  measures,  and  the  30 
resulting  shelf  slope  values  were  assigned  to  the  corresponding 
1-km  interval  sites  within  each  of  the  orientation  segments. 
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Coastal  Process  Variables  (Low  Resolution) 


Wind  Frequency 

The  percentages  of  onshore  winds  and  offshore  winds  greater 
than  two  knots  were  calculated  for  each  of  the  30  orientation 
segments  in  Figure  2,  using  coastal  wind  data  from  stations  at 
Salisbury,  Maryland,  Wallops  Island,  Virginia,  Virginia  Beach, 
Virginia,  Cape  Hatteras,  North  Carolina,  and  Myrtle  Beach,  South 
Carolina  (Garstang  et  al.  1978).  Winds  of  less  than  two  knots 
were  considered  to  have  negligible  sediment  transport  capacity. 
The  data  from  the  nearest  wind  station  were  applied  to  each 
orientation  segment,  and  the  orientation  of  the  best-fit  line 
defining  the  segment  was  used  as  a  criterion  to  determine  which 
wind  directions  to  consider  onshore  and  offshore.  For  each 
segment,  two  of  the  16  compass  directions  were  classed  as 
alongshore,  and  the  other  14  wind  directions  were  divided  evenly 
between  onshore  and  offshore.  Once  the  onshore  and  offshore 
percentages  were  computed  for  a  particular  coastal  segment,  the 
values  for  these  two  variables  were  assigned  to  all  the  1-km 
interval  sites  within  the  orientation  segment. 

Wave  Data 

Wave  data  from  the  Summary  of  Synoptic  Meteorological 
Observations  (S.S.M.O)  were  used  to  calculate  the  percentage  of 
deep-water  waves  greater  than  or  equal  to  1.5-m  high  and  the 
percentage  greater  than  or  equal  to  3.4-m  high  within  each  of 
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the  30  orientation  segments  (U.S.  Naval  Weather  Service  Command, 
1975) .  Offshore-directed  waves  were  not  included  in  the  count. 
Data  from  S.S.M.O.  areas  #15  (Atlantic  City,  New  Jersey)  ,  #16 
(Norfolk,  Virginia),  #17  (Cape  Batteras,  North  Carolina),  and 
#19  (Charleston,  South  Carolina)  were  applied  to  each  orientaton 
segment  within  the  respective  areas.  All  1-km  interval  sites 
within  each  segment  were  given  the  values  for  the  two  wave 
i  variables  that  were  calculated  for  that  segment. 

Tidal  Ranee 

The  mean  tide  range  at  each  study  site  was  estimated  by 
plotting  mean  tide  ranges  at  all  National  Ocean  Survey 
open-coast  tide  stations  (NOAA/National  Ocean  Survey  1982)  on 
1:1, 000, 000-scale  base  maps,  and  interpolating  values  at  all 
sites  located  between  adjacent  tide  stations.  Factors  such  as 

1 

l  inlets,  embayments,  and  major  changes  in  coastline  orientation 

were  taken  into  consideration  in  tide  range  estimations. 

|  IfiJirlsAr-StgrnL-SHrqe,.  flaiaht 

’  Local  values  for  the  storm  surge  heights  with  a  return 

-  period  of  ten  years  were-  extracted  from  three  NOAA  studies  of 

storm  tide  frequency  along  the  mid-Atlantic  coast  (Ho  and  Tracey 

■■ 

’■  1975a.  1975b;  Ho  et  al.  1976).  Storm  tide  height  frequencies 

i 

were  computed  by  NOAA  using:  (1)  the  National  Weather  Service 
|  numerical-dynamic  storm  surge  prediction  model  applied  to  a  full 

1  set  of  climatologically  representative  hurricanes,  and  (2)  tide 


i 
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gage  records  of  winter  (extratropical)  storms  for  locations 


north  of  Cape  Lookout.  North  Carolina.  Hence,  the  ten-year 
storm  surge  height  for  sites  in  Delaware.  Maryland.  Virginia, 
and  North  Carolina  north  of  Cape  Lookout  were  calculated  by 
taking  all  storms  into  consideration.  South  of  Cape  Lookout, 
where  the  relative  frequency  and  magnitude  of  tropical 
cyclone-generated  storm  surges  is  significantly  greater  than 
storm  surges  from  extratropical  storms.  Ho  and  Tracey  (1975a) 
considered  only  hurricanes  in  the  storm-surge  frequency 
analysis.  The  surge  heights  were  computed  in  meters  above  mean 
sea  level. 


Tropical  .Cyglone  .■flfld.Haixigane.  Frequency 
The  number  of  tropical  cyclones  (i.e.r  all  tropical  storms 
and  hurricanes)  and  the  number  of  hurricanes  making  landfall 
between  1886  and  1982  within  92.7-km  (50  nautical-mile)  segments 
of  the  study  area  were  obtained  by  updating  the  totals  reported 
by  Simpson  and  Lawrence  (1971).  Tropical  cyclone  information 
for  1972-1982  was  gathered  from  Neumann  et  al.  (1978)  and  the 
National  Weather  Service's  annual  Atlantic-Car ibbean-Gulf  of 
Mexico  Hurricane  Track  Charts.  The  updating  was  done  by  the 
Simpson  and  Lawrence  (1971)  method,  in  which  tropical  storms 
(sustained  surface  winds  62  to  118  km/hr)  are  counted  only  in 
the  92.7-km  coastal  segment  in  which  the  storm  makes  landfall, 
but  hurricanes  (winds  greater  than  118  km/hr)  are  counted  in  the 
segment  in  which  the  hurricane  makes  landfall  and  in  the 
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adjacent  segment  to  the  right  (i.e. ,  north)  of  landfall.  This 
procedure  takes  into  account  a  hurricane's  destructive  effects 
in  the  right-forward  quadrant  of  its  path  as  the  hurricane  moves 
onshore.  The  tropical  cyclone  frequency  and  hurricane  frequency 
values  obtained  for  each  92.7-km  coastal  segment  were  assigned 
to  all  1-km  interval  sites  within  the  segment. 


REGIONAL-SCALE  SPATIAL  VARIATION 


Initial  analyses  of  the  entire  27-variable  data  set  along 
the  800  km  mid-Atlantic  study  reach  (Table  2)  showed  that  the 
data  matrix  could  be  substantially  reduced  without  significant 
loss  of  information.  This  reduction  was  possible  because  many 
of  the  parameters  have  very  poor  spatial  resolution  (i.e.r 
tropical  cyclone  frequency,  hurricane  frequency,  and  wind 
data) .  Other  variables  provided  unnecessary  resolution  scales 
that  could  be  represented  by  one  inclusive  value  without 
significant  information  loss,  i.e.,  frequency  of  dunes  of 
various  elevations  can  be  represented  by  the  frequency  of  dunes 
above  3  m.  The  3  m  elevation  is  an  approximate  boundary 
between  high  profile  and  low  profile  barrier  islands  in  the 
study  area.  Other  variables  are  merely  standard  deviations  of 
one  of  the  mean  parameters.  The  reduced  data  matrix  includes 
15  variables  (see  Table  4).  Figures  3,  4,  and  5  illustrate  the 
spatial  variation  in  the  15  study  parameters  along  the 
mid-Atlantic  coast  from  Cape  Henlopen,  Delaware,  to  the  North 
Carolina-South  Carolina  border.  Each  of  the  figures  contains 
data  for  five  parameters  and  has  been  divided  into  six 
geographic  sub-regions  corresponding  to  later  discussions  of 
regional  classification.  For  graphical  purposes,  the  data  have 
been  smoothed  by  grouping  the  values  into  classes  representing 
a  range  of  values.  The  actual  data  can  be  found  in  Appendix  A 
for  all  15  variables. 
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Mean  Bar  Number  Mean  number  of  bars  observed  on  +  greater  number  of  bars 

available  air  photos  ~  fewer  bars 
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A  wide  range  of  coastal  orientations  occur  in  the  study 
reach  (Figure  2).  The  reach  between  Cape  Henlopen  and  Ocean 
City  Inlet  is  oriented  dominantly  north-south  with  a  slight 
oceanward  convexity.  From  Ocean  City  Inlet  to  the  Chesapeake 
Bay  the  average  orientation  is  north-northeast,  but  there  are 
departures  from  this  trend  in  the  vicinity  of  inlets, 
particularly  along  the  Virginia  barriers.  The  reach  between 
the  Chesapeake  Bay  and  Cape  Hatteras  strikes  dominantly 
north-northwest,  with  an  oceanward  convexity  occurring  near 
Rodanthe,  North  Carolina,  where  the  strike  becomes  more 
northerly.  The  arcuate  reach  between  Cape  Hatteras  and  Cape 
Lookout  is  oriented  dominantly  northeast,  while  the  arcuate 
reach  between  Cape  Lookout  and  Cape  Fear  strikes  dominantly 
east-northeast.  The  remainder  of  the  southern  North  Carolina 
coastline  strikes  generally  east-west. 

Dune  Frequency 

The  frequency  of  dunes  greater  than  3  meters  in  elevation 
is  an  index  of  barrier  island  topography.  Low-profile  islands 
lack  large  dunes  and  are  subject  to  island-wide  inundation  and 
modification  during  major  overwash  events.  High-profile 
barriers  have  well-established  dune  lines  (largely  man-made) 
and  are  rarely  overwashed  to  the  lagoon  side.  Overwash  on 
high-profile  islands  is  limited  to  breaches  in  the  dune  line 
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and  the  formation  of  small  inter-dune  washover  fans.  Low 
relief  islands  are  subjected  to  periodic  overwash  over  at  least 
one  half  of  the  island  width. 

Between  Cape  Benlopen  and  the  Maryland-Virginia  border, 
high-profile  barriers  are  dominant  with  the  exception  of  small 
regions  at  Cape  Henlopen  and  just  south  of  Ocean  City  along 
northern  Assateague  Island  (Figure  4).  Between  the 
Maryland-Virginia  line  and  the  Chesapeake  Bay  low-profile 
barriers  characterize  the  Virginia  coast.  Fewer  than  25%  of 
the  dunes  along  this  reach  are  greater  than  3  meters  in 
elevation.  The  reach  between  the  Chesapeake  Bay  and  Cape 
Hatteras  is  dominantly  a  high-profile  barrier  coast  with  more 
than  75%  of  the  dunes  above  3  meters  in  most  of  the  area.  The 
only  low-profile  areas  in  this  reach  are  just  south  of  the 
Virginia  border,  near  Oregon  Inlet,  and  midway  between  Rodanthe 
and  Cape  Hatteras.  Almost  the  entire  coastal  reach  between 
Cape  Hatteras  and  Cape  Lookout  is  dominated  by  low-profile 
barriers  with  less  than  25%  of  the  dunes  exceeding  3  meters  in 
elevation.  Between  Cape  Lookout  and  Cape  Fear  the  coast  is 
dominated  by  high-profile  barrier  islands  where  greater  than 
75%  of  the  dunes  are  in  excess  of  3  meters  in  elevation.  The 
southernmost  reach  of  the  study  area  between  Cape  Fear  and  the 
North  Carolina-South  Carolina  border  exhibits  a  highly  variable 
topography,  but  tends  toward  a  dominance  of  low-profile 
barriers. 


The  spatial  frequency  of  inlets  is  highly  variable  along 
the  mid-Atlantic  coast  between  Cape  Benlopen  and  the  North 
Carolina-South  Carolina  border  (Figure  4).  The  northern  reach 
between  Cape  Henlopen  and  Chincoteague  Inlet  and  the  central 
reach  between  the  Chesapeake  Bay  and  Cape  Batteras  contain  few 
inlets.  The  nearly  20  0  km  of  the  coastline  between  Chesapeake 
Bay  and  Cape  Batteras  has  only  one  inlet.  Oregon  Inlet.  Two 
segments  or  the  coast  are  heavily  dissected  by  inlets:  1)  the 
Virginia  barriers  between  Chincoteague  Inlet  and  the  Chesapeake 
Bay,  and  2)  the  southernmost  reach  from  New  River  Inlet  to  the 
South  Carolina  border.  The  reach  between  Cape  Batteras  and  New 
River  Inlet  has  an  intermediate  frequency  of  inlets. 


OPDX  is  not  always  a  true  measure  of  overwash  penetration 
distance.  For  example,  on  high-relief  barriers  overwash 
penetration  is  limited  by  the  dune  barricade.  In  some  cases  on 
low-relief  barriers,  the  entire  barrier  island  is  overwashed, 
hence,  the  OPDX  value  becomes  equal  to,  and  is  dependent  upon, 
island  width.  OPDX  is  essentially  a  measure  of  distance  from 
the  shoreline  to  the  demarcation  between  active  sand  and 
well-established  vegetation.  Therefore,  OPDX  provides  a  good 
measure  of  overwash  penetration  distance  for  the  last  major 
overwash  event  on  low-profile  barriers  that  were  not  totally 
overwashed  to  the  lagoon  (provided  there  has  not  been  total 


revegetation  since  the  event)  .  For  low-profile  barriers  that 
are  totally  overwashed,  OPDX  is  a  conservative  estimate  of 
overwash  distance.  For  high-profile  barriers,  overwash  occurs 
only  at  inter-dune  breaches.  OPDX  is  mapped  in  Figure  3. 

Rate  of  Shoreline  Change  (RSLXl 

RSLX,  the  rate  of  shoreline  change,  is  a  sensitive  measure 
of  the  dynamic  sediment  balance  along  the  coast  during  the 
period  of  aerial  photographic  coverage.  The  majority  of  the 
mid-Atlantic  barriers  are  experiencing  net  shoreline  erosion 
(Figure  3).  The  highest  rates  of  shoreline  erosion  occur  along 
the  Virginia  barriers  between  Chincoteague  Inlet  and  the 
Chesapeake  Bay.  Other  local  areas  of  high  erosion  include:  1) 
the  northern  fifth  of  Assateague  Island  (probably  caused  by 
blockage  of  southerly-drifting  sand  by  engineering  structures 
■  at  Ocean  City  Inlet?  2)  several  zones  of  the  North  Carolina 

coast  just  south  of  the  Virginia-North  Carolina  border?  3)  the 
area  near  Oregon  Inlet?  4)  the  area  just  north  of  the  Carolina 

,  Capes:  Hatteras,  Lookout,  and  Fear?  and  5)  very  localized 

] 

i  areas  on  the  up-drift  sides  of  inlets.  Areas  of  marked  net 

i  accretion  are  few  along  the  mid-Atlantic  barriers, 

i  Rapidly-accreting  areas  are  localized  down-drift  of  inlet 

i  margins  and  on  the  down-littoral  drift  sides  of  some  of  the 

|  capes. 
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Tidal  range  along  the  mid-Atlantic  coast  varies  between 
0.9  m  and  1.6  m,  therefore,  the  entire  study  area  is  classified 


as  microtidal  (Davies  1964).  These  tidal  ranges  are 
significantly  higher  than  Gulf  Coast  tidal  ranges  which  are 
typically  less  than  1  m.  Tidal  range  is  minimal  near  the 
center  of  the  800  km  study  reach  (Figure  4),  and  fluctuates 
between  0.9  m  and  1.1  m  from  the  Chesapeake  Bay  to  New  River 
Inlet.  Tidal  range  increases  gradually  to  the  north  and  south 
from  the  central  region.  To  the  north,  tidal  range  peaks  at 
1.2  m  and  to  the  south  it  increases  to  1.6  m  at  the  North 
Carolina/South  Carolina  border. 

Storm  Surge 

The  storm  surge  values  represent  the  maximum  water  levels 
produced  by  cyclonic  activity  with  a  recurrence  interval  cr  ten 
years.  Storm  surge  in  the  northern  half  of  the  study  area  is 
caused  principally  by  extratropical  storms,  while  southern 
storm  surge  reflects  the  input  of  tropical  cyclones  and 
hurricanes.  The  pattern  of  storm  surge  variation  closely 
correlates  with  the  pattern  of  tidal  range  (Figure  4).  With 
the  exception  of  a  25-km  reach  bordering  the  Chesapeake  Bay  the 
storm  surge  is  lowest  along  the  central  two-thirds  of  the  study 
area  and  increases  to  the  north  and  south.  The  maximum 
northern  storm  surge  values  reach  2  m  along  the  Delaware  coast. 


South  of  Cape  Lookout  the  storm  surge  values  rise  to  2.3  m  near 
Cape  Fear. 
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Mean  grain  size  of  mid-Atlantic  barrier  beaches  ranges 
from  fine  sand  (2.754>,  0.15  mm)  to  coarse  sand  (0.254>,  0.84 
mm) .  Sediment  size  is  highly  variable  throughout  the  reach 
(Figure  4).  This  probably  reflects  a  high  degree  of  dependence 
on  sediment  heredity  and  local  environmental  hydrodynamics. 
Few  spatially  significant  trends  appear  in  the  grain  size  data. 
A  large  area  of  anomalously  fine  sand  occurs  along  the  Virginia 
barriers  and  in  scattered  local  areas  along  the  reach  south  of 
Cape  Lookout.  Areas  of  anomalously  coarse  sand  occur  near 
Duck,  between  Rodanthe  and  Cape  Hatteras,  and  near  Cape  Fear. 

Offshore  Slope 

Offshore  slope  was  measured  from  the  shoreline  to  the 
5.5-m  water  depth  (3-fathom  contour)  and  to  the  9.1-m  depth 
(5-fathom)  contour.  These  measurements  indicate  that  there  is 
a  large-scale  spatial  alternation  of  steep  and  gentle  offshore 
slopes  on  a  scale  of  100  to  200  km  along  the  mid-Atlantic  coast 
(Figure  5).  Offshore  slope  measured  to  the  5.5-m  water  depth 
varies  from  0.93  m/km  to  137.16  m/km  in  the  study  region.  The 
reach  between  Cape  Henlopen  and  the  Maryland-Virginia  border 
has  relatively  steep  offshore  slopes  mostly  in  the  range  of  30 
m/km  to  40  m/km.  Toward  the  southern  end  of  this  reach  the 


values  are  in  the  20  m/km  to  30  m/km  range.  The  Virginia 
barriers  exhibit  the  most  gentle  offshore  slopes  in  the  study 
area,  with  slope  values  mostly  in  the  range  of  2  m/km  to  6 
m/km.  The  reach  between  the  Chesapeake  Bay  and  Cape  Lookout  is 
characterized  by  moderately  gentle  offshore  slopes.  There  is .a 
considerable  degree  of  variation  along  this  reach  but  most  of 
the  slope  values  are  in  the  11  m/km  to  20  m/km  range.  South  of 
Cape  Lookout  offshore  slopes  are  moderately  steep  with  values 
mostly  in  the  range  of  15  m/km  to  30  m/km.  Throughout  the 
study  area  offshore  slopes  proximal  to  inlets  are  lower  due  to 
the  effect  of  offshore  ebb  delta  platforms. 

A  fairly  good  correlation  exists  between  offshore  slope, 
measured  to  the  9.1-m  water  depth,  with  that  measured  to  the 
5.5-m  depth.  The  major  deviations  occur  in  magnitude  of  the 
slope  variations  along  the  coast.  The  Cape  Henlopen  to 
Maryland-Virginia  line  reach  exhibits  only  moderately  steep 
slopes  to  the  9.1-m  contour,  mostly  in  the  range  of  10  m/km  to 
20  m/km.  Larger  areas  of  the  Chesapeake  Bay  to  Cape  Lookout 
reach  are  classified  as  having  gentle  offshore  slopes  with 
values  between  5  m/km  and  10  m/km  when  measured  to  the  9.1-m 
water  depth. 

Island  width 

The  width  of  mid-Atlantic  barrier  islands  ranges  from  less 
than  0.5  km  to  5  km  with  average  widths  between  1  km  and  2  km 
(Figure  3).  Island  width  is  one  of  the  most  variable 
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parameters  observed  along  the  coast  and  is  affected  by  many 
factors  such  as:  1)  proximity  to  inlets;  2)  inheritance  of 
earlier  Holocene  attached  islands;  3)  offshore  slope;  4) 
overwash  processes;  and  5)  inlet  history.  Islands  north  of 
Cape  Hatter as  are  generally  wider  than  those  to  the  south.  The 
widest  islands  are  Assateague,  several  of  the  Virginia 
barriers,  and  the  islands  between  the  Virginia-North  Carolina 
border  and  Oregon  Inlet. 

Lagoon  Width 

Lagoon  width  is  also  highly  variable  like  island  width  and 
ranges  from  zero  to  48.5  km  (Figure  3).  Areas  lacking  lagoons 
occur  along  two  reaches  of  the  Delaware  coast,  along  the 
southern  Virginia  coast,  and  in  localized  areas  of  southern 
North  Carolina.  The  Outer  Banks  between  Oregon  Inlet  and  Core 
Banks  exhibit  the  widest  lagoons  exceeding  15  km.  Areas  of 
moderately  wide  lagoons  (5  km  to  15  km)  occur  along  Assateague 
Island  and  along  the  North  Carolina  coast  between  the 
Virginia-North  Carolina  border  and  Oregon  Inlet.  Narrow 
lagoons  (0.01  km  to  5  km)  occur  along  the  Maryland-Delaware 
coast  north  of  Ocean  City,  in  southern  Virginia,  along  southern 
Core  Banks,  and  between  Cape  Lookout  and  the  South  Carolina 
border. 

Wave  Frequency 


The  frequency  of  waves  greater  than  1.5  m  and  greater  than 


3.4  m  are  highly  correlated  (Figure  5).  In  general,  the 
central  part  of  the  mid-Atlantic  reach  between  the 
Haryland-Virginia  border  and  Oregon  Inlet  has  the  lowest 
frequency  of  waves  of  both  sizes.  Wave  frequency  increases 
northward  and  southward,  with  a  slight  decrease  in  wave 
frequency  over  the  southernmost  70  km  of  the  study  reach. 

Mean  Bar  Number 

Mean  offshore  bar  number,  based  on  aerial  photograph 
observations,  is  highly  variable  in  a  spatial  and  temporal 
sense  along  the  mid-Atlantic  coast  (Figure  5) .  General  trends 
are  apparent  in  bar  configurations.  Most  of  the  Delaware  coast 
is  devoid  of  offshore  bars  in  five  sets  of  photographs 
examined.  Most  of  the  Maryland  coast  oscillates  temporally 
between  no  bars  and  one  bar.  Over  50%  of  the  Virginia 
coastline  between  the  Maryland  border  and  the  Chesapeake  Bay 
has  a  stable  one-bar  system,  but  there  is  considerable  spatial 
variation  between  no  bars  and  one  bar  along  the  Virginia 
barrier  islands  section.  The  reach  with  the  most  spatial  and 
temporal  variability  in  bar  number  is  between  the  Chesapeake 
Bay  and  Rodanthe  (Figure  5) .  Mean  bar  number  varies  between  no 
bars  and  two  bars  in  this  area;  few  sections  are  temporally 
consistent  between  aerial  photo  sets.  A  nearly  stable  one-bar 
pattern  exists  in  all  photo  sets  examined  from  Rodanthe  to  the 
South  Carolina  border. 

Bar  morphologies  were  examined  by  observing  attachment 
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styles.  Figure  6  illustrates  the  range  of  bar  morphologies 
observed  along  the  mid-Atlantic  coast.  The  northern  half  of 
the  Maryland  coast  is  characterized  by  bar  patterns  that 
oscillate  temporally  between  no  bars,  attached  bars,  and  free 
(shore-parallel)  bars.  The  southern  two-thirds  of  Assateague 
Island  is  characterized  by  a  free-bar  system.  Tremendous 
temporal  and  spatial  variability  in  attachment  styles 
characterize  the  Virginia  barriers.  The  reach  between  the 
Chesapeake  Bay  and  Duck,  North  Carolina,  is  characterized  by  a 
stable  free-bar  system.  From  Duck  to  Oregon  Inlet  and  Cape 
Lookout  a  free-bar  system  predominates.  The  northern  half  of 
the  arc  between  Cape  Lookout  and  Cape  Fear  is  characterized  by 
attached  bars,  while  the  southern  half  has  a  free-bar  system. 
Between  Cape  Fear  and  the  North  Carolina-South  Carolina  border 
an  ittached-bar  system  predominates. 
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Figure  . 6.  Morphology  of  longshore  bars 
along  the  mid- Atlantic  coast. 
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DATA  ANALYSIS 


Table  4  summarizes  the  15  parameters  analyzed  along  the 
study  reach,  their  definitions,  and  interpretations  of  negative 
and  positive  variations.  Spatial  relationships  between  the  15 
parameters  were  recognized  using  two  analytical  techniques. 
First,  the  data  set  was  analyzed  by  linear  regression  techniques 
using  the  Statistical  Package  for  the  Social  Sciences  program 
(SPSS)  (Nie  et  al.  1975).  The  results  of  the  regression  analyses 
were  combined  with  our  understanding  of  process  to  establish  a 
hierarchy  of  spatial  relationships  and  trends  visible  along  the 
study  reach.  The  second  mode  of  analysis  was  Principal 
Components  Analysis  (PCA) ,  PCA  was  used  to  look  for  spatial 
organization  and  variations  between  the  15  parameters  along  the 
study  reach.  A  summary  of  the  results  of  each  analysis  is 
presented  in  the  following  sections. 

During  the  initial  analyses  it  became  apparent  that  there 
were  varying  degrees  of  organization  in  the  parameters  depending 
on  the  spatial  scales  being  observed.  Therefore,  analyses  were 
run  for  the  entire  800  km  reach  and  also  for  selected  regional 
data  subsets.  The  study  area  was  divided  into  six  geographic 
subregions  (Figure  2)  according  to  geomorphic  controls  as 
follows:  1)  Cape  Henlopen  to  Chincoteague  Inlet  (115  km  from 
transect  1149-23  to  1116-11);  2)  the  Virginia  barriers  between 
Chincoteague  Inlet  and  the  Chesapeake  Bay  (111  km  from  transect 
1116-11  to  1078-13);  3)  Chesapeake  Bay  to  Cape  Ratteras  (195  km 
from  transect  1077-19  to  1020-21);  4)  Cape  Batteras  to  Cape 
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Lookout  (122  km  from  transect  1019-34  to  983-21);  5)  Cape  Lookout 

to  Cape  Fear  (190  km  from  transect  983-11  to  926-27);  and  6)  Cape 

Fear  to  the  North  Carolina-South  Carolina  border  (56  km  from 

926-17  to  910-10).  Base  data  entered  into  the  analyses  were 

originally  on  a  1  km  spacing  (see  Appendix  A). 

Inlets  and  capes  may  notably  skew  data  for  many  of  the 

classification  parameters  in  sample  sites  peripheral  to  these 

features.  Our  earlier  studies  of  spatial  variation  in  rates  of 

shoreline  change  suggest  that  the  along-the-coast  extent  of  cape 

influence  is  about  4  km  and  the  along-the-coast  extent  of  inlet 

influence  is  less  than  2  km  (Vincent  et  al.  1976;  Dolan  et  al. 

1977).  Therefore,  each  regression  and  principal  components 

analysis  was  run  for  data  sets  including  and  excluding  inlet  and 
• 

cape  effects.  Data  sets  excluding  inlet  and  cape  effects 
excluded  samples  within  2  km  of  inlets  and  within  4  km  of  capes. 

CoixelaUan.  .Analysis 

This  section  contains  summaries  of  the  significant 
correlations  between  the  15  classification  parameters  described 
in  Table  4.  Only  correlations  significant  at  the  level  of  * 
0.001  (99.99%)  were  accepted.  Significant  values  of  r  were 
determined  using  the  Test  Statistic  (T) t 
r  /  n  -  2 

T  -  _ 

/  l  -  rZ 

Values  for  T  were  taken  from  the  table  in  Kleinbaum  and  Kupper 
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(1976).  Values  of  r  ■  0.35  were  chosen  as  lower  limits  of 
acceptance  as  long  as  they  were  above  the  minimum  calculated 
level  of  significance.  Discussions  of  correlations  are 
stratified  as  moderate  (r  between  0.35  and  0.59)  and  strong  (r 
greater  than  0.60).  Only  three  data  sets  contained  few  enough 
cases  to  mandate  using  0.60  as  the  minimally  significant  r  value 
(see  data  sets  REGBRM2 ,  REGFWI,  and  REGFRM2 )  in  Table  6. 

Correlation  analyses  were  used  to  answer  a  number  of 
questions  including:  1)  What  associations  exist  between 
variables  and  what  dependencies  are  suggested  by  these 
associations?  2)  Are  the  associations  observed  for  the  entire 
coast  also  observed  when  selected  subregions  are  studied,  i.e., 
are  there  different  patterns  and  scales  of  organization  between 
these  variables  along  the  coast?  3)  Does  the  presence  of  capes 
and  inlets  affect  the  associations  observed  between  variables? 
The  correlation  analysis  does  not  tell  us  what  variables  are 
important  in  terms  of  classifying  the  coast  because  correlation 
analysis  does  not  deal  directly  with  the  questions  of  how  the 
relationships  between  the  coastal  parameters  behave  in  a  spatial 
manner  or  along  the  coast  or  what  their  relative  magnitudes  are 
in  given  regions.  These  latter  questions  will  be  addressed  using 
Principal  Components  Analysis  in  a  later  section  of  this  report. 

Entire  Coast  -  Cape  Henlopen,  Delaware,  to  North  Carolina- 

South  Carolina  Border 

In  this  data  set  800  cases  provided  the  base  data  at  1-km 
intervals  from  Cape  Henlopen,  Delaware,  to  the  North 


Carolina-South  Carolina  border.  Five  different  correlations  were 
run  from  this  data  set  as  follows:  1)  BIGDAT  -  the  entire  coast 
at  1-km  intervals,  including  areas  peripheral  to  inlets  and  capes 


(n  •  800)  ;  2)  INLETR2  -  the  entire  coast  at  1-km  intervals, 
excluding  areas  peripheral  to  inlets  and  capes  (n  -  564);  3) 
SUBSETS  -  the  entire  coast  at  5-km  intervals,  including  areas 
peripheral  to  inlets  and  capes  (n  ■  160) ;  4)  BIGDAT5  -  the  entire 
coast  averaged  over  5-km  sections  (n  *  160) ;  and  5)  INLTRM5  -  the 
entire  coast  averaged  over  5-km  sections,  excluding  areas 
peripheral  to  inlets  and  capes  (  n  ■  118).  Significant  moderate 
and  strong  correlations  are  summarized  in  Table  5. 

The  Entire  Coast  at  1-km  Intervals  (BIGDAT  and 

INLETR2? 

Correlation  analysis  of  the  15  variables  for  the  entire 
coast  at  1-km  intervals  shows  strong  correlations  between  storm 
surge  and  tidal  range  and  between  shoreline  strike  and  the 
frequency  of  waves  greater  than  1.5  m.  Moderate  correlations 
occur  between  the  following  variables:  coastal  strike,  frequency 
of  dunes  higher  than  3  m.  spatial  inlet  frequency,  overwash 
penetration  distance,  tidal  range,  storm  surge,  sediment  size, 
offshore  slope  to  the  5.5-m  depth,  offshore  slope  to  the  9.1-m 
depth,  island  width,  lagoon  width,  frequency  of  waves  greater 
than  1.5-m  high,  frequency  of  waves  greater  than  3.4-m  high,  and 
mean  bar  number  (see  Table  5).  Differences  do  occur  in  the 
numbers  of  significant  correlations  in  data  sets  including  sites 
adjacent  to  inlets  and  capes  compared  with  data  sets  excluding 


TABLE  5 

SUMMARY  OF  SIGNIFICANT  CORRELATIONS  Cn  -  0.001  level) 
CAFE  HENLOPEN  TO  THE  NORTH  CAROLINA-SOUTH  CAROLINA  BORDER 


Parameter4 

BIGDAT 

INLETR2 

BIGDAT5 

SUBSET5 

INLTRM5 

STRK 

WPQK-.50) 

WFQ3 (-.36) 

WPQK-.61)* 

WFQ3 (-.48) 

DFQ3 (+.39) 
WPQM-.52) 

WFQ3 (-.37) 

WPQK-.52) 

WFQ3 (-.37) 

DFQ3 (+.38) 
WFQ1 (-.63) * 
WFQ3 (-.52) 

DFQ3 

OPDX (-.35) 

OPDX(-.39) 

OPDX (-.40) 

OFS5 (+.40) 
OFS9(+.36) 

STRK (+.39) 

OPDX (-.37) 

OPDX  (-.41) 
STRK(+.38) 

INFQ 

SEDS (-.43) 

TDRG(+.37> 

STSG (+.41) 

SEDS (-.38) 

SEDS (-.53) 

OFS9 (-.40) 

SEDS (-.42) 

TDRG(+.39) 
STSG (+.45) 
SEDS(-.44) 

OPDX 

LAGW(+.35) 

DFQ3 (-.35) 

DFQ3 (-.39) 

LAGW(+. 40) 

DFQ3 (-.40) 

DFQ3 (-.37) 

DFQ3 (-.41)  | 

RSLX  - - - - - - - — — - 

— i 

\ 

1 

TDRG 

STSG (+.64)* 

STSG(+.68)  * 

INFQ (+.37) 

STSG (+.65)* 

STSG (+.61)* 

STSG(+.69) 

INFQ(+,39) 

4 

STSG 

LAGW (-.40) 
TDRG(+.64) • 

ISLW(-.35) 

LAGW (-.40) 

INFQ(+ .  41) 

TDRG (+.68)* 

LAGW (-.41) 
TDRG(+.65) * 

LAGW (-.41) 
TDRG(+. 61) * 

ISLW(-.43) 
LAGW  (-.40) 
INFQ (+.45) 
TDRG(+.69) 

* 

SEDS 

INFQ(-. 43) 

INFQ(-.38> 

LAGW (+.35) 

INFQ (-.53) 

INFQ(-.42) 

LAGW (+.35) 
INFQ(-.44) 

OFS5 

WFQ3(+.45) 

WFQ3 (+.47) 
BARS(-.44) 

WFQ3 (+.55) 
BARS(-.37) 

DPQ3 (+.40) 

WFQ3 (+.37) 

ISLW (-.36) 
WFQ3 (+.53) 
BARS (-.53) 

OFS9 

WFQ3 (+.41) 

WFQ3 (+.48) 

WFQ3 (+.53) 

DFQ3 (+.36) 

INFQ (-.40) 

WFCttt+.SS)1 

I 

ISLW 

WFQK-.45) 

WFQ3 (-.36) 

STSG (-.35) 

WPQ1(-.56)J 
WFQ3 (-.44) 
STSG (-.43) 
OFS5(-.36>; 

LAGW 

OPDX (+.35) 
8TSG(-.40) 

STSG(-.40) 

OPDX (+.40) 
STSG(-.41) 

SEDS (+.35) 

STSG (-.41) 

STSG (-.40) 
SEDS (+.35) 

WFQ1 

STRK (-.50) 

STRK (-.61) * 

ISLW (-. 45) 

STRK (-.52) 

STRK (-.52) 

STRK (-.63* 
ISLW  (-.5611 

WFQ3 

8TRK(-.36) 

OF85 (+.45) 
OF89(+.41) 

BARS (-.3 9) 

STRK (-.48) 

OFS5 (+.47) 

OFS9 (+.48) 

ISLW (-.36) 

STRK (-.57) 

OFS5 (+.55) 

OFS9 (+. 53) 

STRR(-.37) 

BARS (-.43 
STRK (-.52 
OFSS (+.53 
OFS9 (+.5S 
ISLW (-.44 

BARS 

OFS5 (-.44) 

BARS (-.39) 

OFSS (-.37) 

OrS5(-.5J 
WFQ3 (-.43 

"  aanotta  particularly  strong 

correlations. 

*  for  parameter  definitions  see  Table  4. 


these  values. 

Several  correlations  were  common  to  both  data  sets  which 
illustrates  their  persistence  along  the  coast  and  independence 
from  the  effects  of  cape  and  inlet  processes.  The  strongest 
correlation  was  between  tidal  range  and  storm  surge.  As  tidal 
range  increases  storm  surge  also  increases.  The  frequency  of 
large  waves  approaching  the  coast  is  related  to  the  strike  of  the 
coastline.  As  the  coast  strikes  in  a  more  easterly  direction, 
the  frequency  of  waves  above  1.5  m  and  above  3.4  m  decreases. 
Offshore  slope  also  appears  to  directly  affect  the  frequency  of 
large  waves  (above  3.4  m)  where  large  waves  occur  more  frequently 
in  areas  of  steeper  offshore  slope.  Overwash  penetration 
distance  is  controlled  by  island  topography.  As  the  frequency  of 
dunes  above  3  meters  increases,  OPDX  decreases. 

Additional  correlations  appear  when  the  sites  adjacent  to 
inlets  and  capes  are  removed,  thereby  filtering  out  the  direct 
effects  of  these  features  (data  set  INLETR2) .  Inlet  frequency 
appears  to  be  related  to  tidal  range  and  storm  surge.  As  the 
frequency  of  inlets  along  the  coast  increases,  tidal  range  and 
storm  surge  increase.  Island  width  decreases  as  storm  surge 
increases.  Mean  bar  number  appears  to  be  related  to  the 
frequency  of  large  waves  and  to  offshore  slope.  The  number  of 
bars  decreases  as  offshore  slope  increases  and  as  the  frequency 
of  waves  greater  than  3.4  m  increases.  These  variables  affecting 
mean  bar  number  are  masked  in  the  vicinity  of  inlets  due  to  the 
complex  local  hydrodynamics  in  these  areas. 
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The  Entire  Coast  at  5-km  Intervals  (BIGDAT5) 

At  the  outset  of  this  study  we  decided  that  a  1-km  sampling 
interval  may  be  required  to  detect  regional  and  local  trends  and 
organizational  patterns  in  the  study  variables.  BIGDAT5 
represents  a  subset  of  the  entire  80  0  case  data  set  composed  of 
every  fifth  case,  hence  this  represents  a  sampling  interval  of  5 
km  along  the  coast.  The  relationships  in  the  medium  resolution 
sampling  scheme  of  5  km  (BIGDAT5)  are  compared  with  high 
resolution  sampling  (1  km)  shown  in  BIGDAT  and  INLETR2  (see  table 
5).  In  this  manner,  the  persistence  of  the  relationships  can  be 
tested  for  .different  sampling  resolutions. 

There  is  not  a  direct  correspondence  between  trends  visible 
in  the  two  data  sets.  Over  75%  of  the  same  significant 
correlations  occur  in  both  (Table  5) ,  but  there  are  numerous 
strong  relationships  that  occur  in  only  one  of  the  data  sets. 
The  following  correlations  occur  solely  in  the  INLETR2  data  set: 
1)  as  inlet  frequency  increases,  tidal  range  and  storm  surge 
increase;  2)  as  storm  surge  increases,  island  width  decreases;  3) 
as  wave  frequency  increases,  island  width  decreases;  and  4)  as 
wave  frequency  above  3.4  m  increases,  the  mean  number  of  bars 
increases.  The  correlations  limited  to  the  BIGDAT5  data  set  are: 
1)  as  coastal  strike  becomes  more  easterly,  the  frequency  of 
dunes  above  3  m  increases;  2)  as  dune  frequency  increases, 
offshore  slope  increases,  3)  inlet  frequency  increases  as 


offshore  slope  to  the  9.1-m  depth  decreases;  4)  as  lagoon  width 
increases ,  OPDX  and  sediment  size  also  increase. 


From  this  data  it  is  unclear  what  the  effects  of  different 
sampling  resolutions  are  on  the  outcome  of  the  correlations.  The 
coarser  sampling  interval  failed  to  pick  up  correlations  with 
wave  frequency  and  storm  surge#  which  are  relatively  low 
resolution  parameters.  Therefore#  a  strong  case  is  made  for  the 
use  of  the  finer  sampling  interval,  and  the  correlations  obtained 
with  the  1-km  interval  analysis  are  probably  more  reliable. 

Entire  Coast  Averaged  Over  5-km  Segments  (SUBSET5 

and  INLTRM57 

Considerable  differences  exist  in  the  numbers  of  significant 
correlations  between  5-km  averages  taken  including  sites  adjacent 
to  inlets  and  capes  (SUBSET5 )  and  5-km  averages  with  inlet  and 
cape  sites  removed  (Table  5) .  Correlations  in  S0BSET5  and  BIGDAT 
are  very  similar.  The  only  differences  are  that  the  BIGDAT  data 
set  shows  a  positive  correlation  between  OPDX  and  lagoon  width 
and  a  positive  correlation  between  the  frequency  of  large  waves 
and  offshore  slope,  while  no  significant  correlations  occurred 
between  these  variables  when  the  entire  coast  is  averaged  over  5 
km.  The  strong  correspondence  between  SUBSET5  and  BIGDAT  and  the 
weaker  correspondence  between  BIGDAT5  arid  BIGDAT  indicates  that 
while  a  coarsening  of  sample  interval  from  1  km  to  5  km  results 

in  a  loss  of  information#  a  smoothing  of  the  1-km  data  by 
averaging  over  5  s.n  intervals  does  not  significantly  affect  the 
results  of  the  correlations.  Likewise#  there  are  very  few 
differences  in  correlations  between  INLETR2  and  INLTRM5#  which 


suggests  that  a  similar  smoothing  of  high  resolution  sample  data 
without  inlet  and  cape  sites  can  be  done  without  disturbing  the 
final  relationships. 


Effects  Q-f-  Inlets  and  Capes  on  Correlations 
Removing  the  cases  adjacent  to  inlets  and  capes  does  affect 
correlations  in  the  1-km  and  5-km  sampling  schemes.  The  major 
effect  of  removing  sites  proximal  to  inlets  and  capes  appears  to 
be  an  increase  in  the  value  of  the  correlation  coefficients 
observed  in  the  original  data  sets.  Of  secondary  importance, 
several  new  associations  appeared  that  were  not  observed  in  the 
original  data  sets  using  all  of  the  cases.  Host  of  these  new 
associations  were  apparent  in  the  original  data  sets  but  their 
correlation  coefficients  were  just  below  the  0.001  level  of 
significance.  Removing  the  cases  adjacent  to  inlets  and  capes 
has  removed  a  significant  amount  of  noise  from  the  system  and 
allowed  the  associations  to  be  more  readily  observed.  The  major 
new  associations  observed  are  the  following:  1)  positive 
correlations  of  inlet  frequency  with  storm  surge  and  tidal  range; 
2)  negative  correlation  between  island  width  and  wave  frequency, 
and  3)  negative  correlation  between  offshore  slope  and  bar 
number. 

Geographic  Subregions  at  1-km  Sampling  Intervals 
Subsets  of  the  800-km  data  set  were  created  in  accordance 
with  major  geomorphic-geographic  boundaries  along  the 
mid-Atlantic  coast  in  order  to  determine  if  there  was  any 
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significant  organization  between  the  various  parameters  on  a 
regional  scale  that  may  be  masked  by  analyses  of  the  entire 
coast.  Table  6  shows  the  significant  correlations  for  the 
various  geographic  subregions.  Table  6  clearly  shows  that 
greater  numbers  of  significant  correlations  as  well  as 
correlations  with  higher  r  values  occur  for  individual  geographic 
subregions  compared  to  the  entire  coast.  This  indicates  that 
many  of  the  relationships  between  parameters  are  organized  on  a 
region-specific  scale  that  changes  its  overall  pattern  along  the 
coast.  In  addition,  there  are  differences  in  the  degree  of 
correlations  between  the  various  subregions.  For  example, 
subregion  A  between  Cape  Henlopen  and  Chincoteague  Inlet  has  the 
highest  number  of  correlations,  while  other  subregions  have  fewer 
significant  correlations.  Therefore,  subregion  A  appears  to  be 
structured  in  a  more  orderly  fashion  with  respect  to  the 
variables  studied.  Subregion  B  shows  the  lowest  number  of 
associations  between  variables,  suggesting  that  this  is  the  least 
structured  subregion  along  the  entire  mid-Atlantic  coast. 

Another  general  trend  observed  in  the  regional  data  sets  is 
that  subregions  with  consistent  coastal  strike  tend  to  exhibit 
higher  degrees  of  organization  than  regions  with  major  shifts  in 
orientation.  This  suggests  that  coastal  strike  is  a  major 
independent  variable  controlling  the  development  of  many  of  the 
other  variables  such  as  wave  frequency,  offshore  slope,  and  storm 
surge.  In  general,  the  major  trends  observed  in  the  entire 
coastal  data  sets  become  much  less  obvious  in  the  subregional 
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data  sets.  For  example,  the  negative  correlation  between  inlet 
frequency  and  sediment  size  that  was  observed  in  all  five  of  the 
entire  coastal  data  set  appears  in  only  two  of  the  six 
subregional  data  sets.  Relationships  between  coastal  strike  and 
wave  frequency  are  reversed  from  what  was  observed  in  the  entire 
coastal  data  set. 

A  large  number  of  new  associations  are  observed  in  the 
subregional  data  sets  that  are  not  present  in  correlations  of  the 
entire  coastal  data  sets.  However,  none  of  these  trends  appear 
in  as  many  regions  as  the  major  trends  first  observed  in  analyses 
of  the  entire  data  set.  hence,  they  appear  to  be  specific  to 
individual  regions.  New  associations  common  to  at  least  three  of 
the  six  subregions  include  the  following:  1)  tidal  range 
increases  as  coastal  strike  becomes  more  easterly;  2)  storm  surge 
increases  as  lagoon  width  decreases;  3)  offshore  slope  increases 
as  lagoon  width  decreases;  4)  as  island  width  decreases,  the 
frequency  of  large  dunes  increases;  and  5)  as  the  frequency  of 
large  waves  increases,  the  number  of  offshore  bars  decreases. 

Subregion.  A_ Cape  Henlopen  to  Chlncoteaque  Inlet 

Table  6  shows  that  a  great  number  of  significant 
correlations  occur  in  the  region  between  Cape  Henlopen  and 
Chincoteague  Inlet.  Very  minimal  differences  occur  between  the 
data  set  with  inlets  and  capes  compared  to  the  one  with  these 
areas  excluded.  The  only  significant  correlations  in  subregion  A 
(REGARM)  after  inlet  and  cape  areas  were  excluded  were:  1) 
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negative  correlation  between  dune  frequency  and  island  width;  2) 
positive  correlation  between  dune  frequency  and  more  easterly 
trending  coastal  strike;  3)  negative  correlation  between  OPDX  and 

i 

offshore  slope  to  9.1-m  depth;  4)  positive  correlation  between 
offshore  slope  to  5.5-ra  depth  and  tidal  range;  5)  negative 
correlation  between  sediment  size  and  mean  bar  number;  and  6) 
positive  correlation  between  tidal  range  and  offshore  slope  to 
the  5.5-m  depth. 

The  strongest  correlations  will  be  summarized  briefly.  The 
moderate  correlations  are  numerous  and  can  be  studied  in  Table  6. 
Coastal  strike  appears  to  be  strongly  related  to  tidal  range, 
storm  surge,  and  mean  bar  number.  As  the  strike  of  the  coast 
becomes  more  easterly,  storm  surge  and  tidal  range  increase  while 
mean  bar  number  tends  to  decrease.  The  frequency  of  dunes 
greater  than  3-m  high  increases  as  the  frequency  of  waves  above 
1.5  m  increases.  Inlet  frequency  is  related  to  lagoon  width  and 
mean  bar  number  such  that  inlet  frequency  increases  as  bar  number 
decreases  and  as  lagoon  width  decreases.  Tide  range  and  storm 
surge  are  highly  correlated  in  a  positive  sense.  Storm  surge 
increases  as  offshore  slope  to  the  5.5-m  depth  increases.  The 
mean  number  of  offshore  bars  appears  to  be  controlled  by  coastal 
strike,  inlet  frequency,  storm  surge,  and  offshore  slope. 
Significant  correlations  also  occur  between  offshore  bars  and 
dune  frequency,  OPDX,  tidal  range,  island  width,  and  lagoon 
width.  Mean  bar  number  increases  as  coastal  strike  becomes  more 
easterly,  as  dune  frequency  decreases,  as  inlet  frequency 
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decreases,  as  OPDX  increases,  as  tidal  range  decreases,  as  storm 
surge  decreases,  as  sediment  size  decreases,  as  offshore  slope 
decreases,  as  island  width  increases,  and  as  lagoon  width 
increases. 


Subregion  B  -  Chincoteaoue  Inlet  to  the  Chesapeake  Bav 
Very  few  associations  occur  between  variables  along  the 
Virginia  barriers.  The  apparent  lack  of  geomorphic  organization 
within  this  subregion  is  most  likely  due  to  the  high  frequency  of 
inlets  and  great  local  variance  in  shoreline  orientation.  The 
only  moderate  correlations  are:  1)  sediment  size  tends  to 
increase  as  tidal  range  decreases  and  as  lagoon  width  decreases; 
and.  2)  lagoon  width  tends  to  increase  as  offshore  slope 
decreases,  as  sediment  size  decreases,  as  the  frequency  of  large 
waves  increases,  and  as  tidal  range  increases. 

Subregion  C  -  Chesapeake  Bayi-to— Cape  Hat.texaa 
The  strong  associations  are  as  follows:  1)  as  coastal 
strike  becomes  more  northerly,  lagoon  width  decreases  and  the 
frequency  of  large  waves  increases;  and  2)  storm  surge  increases 
as  inlet  frequency  increases.  A  large  number  of  moderate 
correlations  occur  along  this  reach  of  the  coast  and  are 
summarized  in  Table  6. 
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Subregion  D  -  Cape  Hatteras  to.Cape  Lockout 


Strong  correlations  occur  between  the  following  variables  in 
subregion  D:  1)  as  coastal  strike  becomes  more  easterly,  the 
frequency  of  large  dunes  increases,  storm  surge  decreases, 
offshore  slope  decreases,  and  lagoon  width  increases;  2)  as  storm 
surge  increases,  the  frequency  of  large  dunes  decreases,  tidal 
range  increases,  lagoon  width  decreases,  offshore  slope 
decreases,  and  coastal  strike  becomes  more  northerly;  and  3)  as 
offshore  slope  increases,  coastal  strike  becomes  more  easterly 
and  storm  surge  decreases.  A  large  number  of  significant 
moderate  correlations  occur  along  this  reach  of  the  coast  and  are 
summarized  in  Table  6. 

SUbregion  £  -  Cape  Lookout  to  Cape  Fear 
Strong  correlations  occur  between  the  following  variables  in 
subregion  E:  1)  as  coastal  strike  becomes  more  easterly,  storm 
surge  decreases  and  tidal  range  also  decreases;  2)  as  the 
frequency  of  large  dunes  increases,  the  overwash  penetration 
distance  decreases;  3)  overwash  penetration  distance  increases  as 
island  width  increases,  as  the  frequency  of  large  dunes 
decreases,  and  as  the  number  of  offshore  bars  decreases;  4)  tidal 
range  increases  as  storm  surge  increases;  and  5)  the  frequency  of 
large  waves  increases  as  lagoon  width  decreases.  A  large  number 
of  significant  moderate  correlations  occur  along  this  reach  of 
the  coast  and  are  summarized  in  Table  6. 
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Subregion  F.  -..Cape  Fear  to  the. North  Carolina- 
South  Carolina  Border 

In  region  F  only  moderate  correlations  occur  at  the  0.001 
level  of  significance  due  to  the  small  data  set  after  areas 
adjacent  to  inlets  and  capes  were  removed.  The  moderate 
correlations  are:  1)  as  coastal  strike  becomes  more  easterly, 
inlet  frequency  decreases,  overwash  penetration  distance 
decreases,  tidal  range  decreases,  storm  surge  increases,  offshore 
slope  decreases,  and  island  width  increases;  2)  as  inlet 
frequency  increases,  tidal  range  increases,  storm  surge 
decreases,  and  the  number  of  offshore  bars  decreases;  3)  overwash 
penetration  distance  increases  as  tidal  range  increases,  as 
island  width  increases,  and  as  the  coastal  strike  becomes  more 
northerly;  4)  storm  surge  increases  as  offshore  slope  decreases, 
as  island  width  decreases,  as  coastal  strike  becomes  more 
easterly,  as  inlet  frequency  decreases,  as  overwash  penetration 
distance  decreases,  and  as  tidal  range  decreases;  5)  offshore 
slope  increases  as  lagoon  width  decreases,  as  coastal  strike 
becomes  more  northerly,  and  as  storm  surge  decreases;  and  6)  the 
number  of  bars  increases  as  inlet  frequency  decreases  and  as 
tidal  range  decreases. 
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Principal  components  analysis  is  an  analytical  method 
applicable  to  large  data  matrices  which  transforms  a  series  of 
correlated  variables  into  a  new  set  of  statistically  independent 
(orthogonal)  factors  called  principal  components  or  eigenvectors 
(Kleinbaum  and  Kupper  1978).  The  first  eigenvector  explains  the 
largest  amount  of  variance  in  the  system  while  subsequent 
eigenvectors  explain  successively  smaller  amounts  of  the  total 
variance.  PCA  also  transforms  original  data  scores  into 
weightings  (scores)  where  one  unique  set  of  scores  occurs  with 
each  principal  component  (Daultrey  1976)  .  This  technique  has 
been  successful  in  explaining  the  variance  in  coastal 
geomorphology  and  beach  systems  (Vincent  et  al.  1975;  Winant  and 
Aubrey  1976;  Resio  et  al.  1977,  and  Fisher  et  al.  1982)  and  is 
appropriate  for  the  analysis  of  the  15  variables  in  this  study 
listed  in  Table  4. 

Principal  components  analysis  was  run  for  the  entire  data 
set  (800  cases)  at  1-km  intervals  and  also  for  the  various 
geographic  subregions  denoted  in  Figure  2.  Only  the  runs  where 

cases  adjacent  to  inlets  and  capes  were  excluded  will  be 

\ 

discussed  in  this  report. 

Entire  Coast  from  Cape  Henlopen  to  the  South  Carolina*- 

North  Carolina  Border 

The  data  set  excluding  areas  near  capes  and  inlets  includes 
564  cases.  The  first  four  eigenvectors  are  statistically 
significant  (according  to  Overland  and  Preisendorf er  1982). 
These  four  account  for  631  of  the  total  variance  in  the  data. 


Eigenvector  1  accounts  for  23%  of  the  variance  alone.  Figure  7 
schematically  depicts  the  results  of  the  PC  A  interpreted  with 
respect  to  the  coastal  geomorphology  and  processes.  The 
schematic  models  shown  in  Figure  7  were  constructed  in  the 
following  manner.  First,  the  basic  model  was  developed  using  the 
weightings  of  the  statistically  significant  eigenvectors. 
Second,  the  significant  eigenvectors  were  reconstructed  (merged) 
using  a  program  that  uses  the  mean,  standard  deviation,  and 
weighting  of  each  of  the  eigenvectors  in  accord  with  the 
following: 

Rij  -  Xi(±c)  <ai)  (ei;j) , 

where  R  is  the  reconstructed  value  for  the  ith  variable  and  the 

« 

jth  vector  is  the  mean  -for  each  variable,  c  is  a  constant,  ai 
is  the  standard  deviation  of  each  variable,  and  c^j  is  the 
loading  on  each  vector  for  each  variable.  For  a  detailed 
discussion  of  this  technique  see  Resio  et  al.  (1974).  Analysis 
of  the  entire  coastal  data  set  is  primarily  sensitive  to 
large-scale  regional  variations,  while  analysis  of  regional 
subsets  is  more  likely  to  be  depicting  sroaller-scale  variations 
and  organization  in  the  data.  Discussion  of  the  results  of  PCA 
will  proceed  from  north  to  south  along  the  coast. 

The  northern  reach  of  the  mid-Atlantic  barrier  coast  from 
Cape  Benlopen  to  the  southern  part  of  Assateague  Island, 
Virginia,  is  characterized  by  high-profile  barriers  (with  the 
exception  of  northernmost  Assatague  Island)  with  very  steep 
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offshore  slopes  (the  steepest  slopes  of  the  mid-Atlantic  coast) . 
Few  inlets  occur  along  this  north  to  north-northeasterly  striking 
coastline.  Islands  are  variable  in  width  but  generally  increase 
in  width  toward  the  south.  Lagoon  width  also  tends  to  increase 
from  no  lagoon  in  the  north  to  wide  lagoons  in  the  south.  The 
frequency  of  large  waves  is  high  along  this  reach  while  tidal 
range  and  storm  surge  are  near  the  mean  values  for  the 
mid-Atlantic  coast.  Sediments  are  coarse  to  moderately  coarse 
and  single  bars  predominate.  Overwash  penetration  distance  is 
about  average  and  most  of  the  region  is  experiencing  slight  net 
erosion. 

Islands  between  southern  Assateague  Island  and  the 
Chesapeake  Bay  (Virginia  barriers)  have  low  profiles  and  very 
gentle  offshore  slopes.  These  islands  strike  dominantly 
north-northeast  and  are  dissected  by  numerous  inlets  (the  highest 
inlet  frequency  along  the  mid-Atlantic  coast)  .  The  islands  in 
the  northern  half  of  this  reach  are  relatively  wide  while  the 
southern  islands  are  narrow.  Lagoon  width  generally  increases 
toward  the  south.  The  frequency  of  large  waves  is  very  low  and 
tidal  range  and  storm  surge  are  near  the  regional  mean  for  the 
entire  coast.  Sediments  are  moderately  fine  to  fine  and  a 
single-bar  system  predominates.  Overwash  penetration  distance  is 
high  and  erosion  rates  are  among  the  highest  of  the  entire  coast. 

South  from  the  Chesapeake  Bay  to  the  Kitty  Hawk  area  of 
North  Carolina  (south  of  the  CERC  Pier),  the  barrier  islands  are 
generally  high  profile.  Offshore  slopes  along  this  reach  are 
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moderately  gentle  in  the  north  and  become  steeper  southward.  The 
strike  of  these  islands  is  north-northwest  throughout  this  reach 
and  the  area  has  no  inlets.  Islands  are  wide  along  this  reach 
while  lagoons  are  of  average  width.  The  frequency  of  large  waves 
is  very  low.  Storm  surge  is  high  along  the  northern  half  of  the 
reach  and  low  to  the  south.  Tidal  range  is  low  at  the  north  and 
increases  to  about  the  mean  toward  the  south.  Beach  sediment 
size  is  variable  along  this  reach  but  is  generally  coarse  to  the 
north  and  south  and  fine  along  the  central  portion  of  this  reach. 
A  two-bar  system  predominates  offshore  in  most  areas.  Overwash 
penetration  distance  is  low  along  the  northern  and  southern 
thirds  of  this  reach  and  high  in  the  central  portion  of  this 
slowly  eroding  reach  of  the  coast'. 

Between  the  Kitty  Hawk  area  and  Cape  Hatteras  high-profile 
barriers  are  dominant  except  for  the  area  bordering  Oregon  Inlet. 
Between  Kitty  Hawk  and  Rodanthe,  the  North  Carolina  coast  strikes 
north-northwest.  South  of  Rodanthe  the  coast  strikes  more 
northerly.  Offshore  slope  along  this  reach  is  moderately  steep 
to  steep,  except  for  a  gentle  reach  near  Oregon  Inlet.  Islands 
and  lagoons  are  wide  along  the  entire  reach  and  only  one  inlet 
(Oregon  Inlet)  occurs  in  the  area.  The  frequency  of  large  waves 
is  moderate  near  Oregon  Inlet  and  increases  toward  Cape  Hatteras. 
Storm  surge  values  are  near  the  mean  and  tidal  range  is  low  along 
this  reach.  Beach  sediments  are  moderately  coarse  near  Oregon 
Inlet  and  coarser  toward  Cape  Hatteras.  The  dominant  bar 
morphology  is  a  two-bar  pattern  with  small  areas  of  alternating 
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one-bar  and  two-bar  systems.  Overwash  penetration  distance  is 
moderate  to  great  along  this  reach  and  most  of  this  coastline  is 
experiencing  moderately  rapid  erosion. 

The  region  between  Cape  Batteras  and  Cape  Lookout  is 
characterized  by  north-northeasterly  to  northeasterly  striking 
islands  of  low  profile.  Offshore  slope  is  moderately  steep  to 
moderate.  The  islands  are  wide  along  the  northern  half  of  the 
region  and  narrow  to  the  south,  and  lagoons  are  wide  throughout. 
Few  inlets  occur  in  this  region,  hence  it  has  an  inlet  frequency 
near  the  mean.  The  frequency  of  large  waves  is  high  along  this 
reach,  storm  surge  is  average,  and  tidal  range  increases  from  low 
to  high  southward  along  this  reach.  Sediments  are  variable  along 
this  reach,  but  are  near  the  mean  size  for  the  mid-Atlantic 
barrier  beaches.  The  northern  half  of  the  area  is  dominated  by 
two  bars  while  the  southern  half  is  dominated  by  alternating 
one-bar  and  two-bar  systems.  Overwash  penetration  distance  is 
high  along  the  entire  area  and  rates  of  shoreline  change  are  low. 

From  Cape  Lookout  south  toward  Cape  Fear  coastal  strike 
progressively  shifts  from  east-southeast  to  north-northeast. 
Most  of  the  reach  is  dominated  by  high-profile  barriers  except 
for  the  northern  segment  between  Beaufort  Inlet  and  Cape  Lookout 
and  the  southernmost  segment  between  Fort  Fisher  and  Cape  Fear. 
This  reach  of  the  coast  has  a  moderately  steep  offshore  slope 
except  for  a  gently-sloping  area  between  Bogue  Inlet  and  central 
Ashe  Island.  The  islands  are  relatively  narrow  throughout  this 
reach  as  are  the  lagoons,  except  for  an  area  of  wider  lagoons 
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north  of  Bogue  Inlet.  Few  inlets  occur  along  the  northern  half 
of  this  reach  while  many  inlets  occur  along  the  southern  half. 
Large  waves  are  very  frequent  along  roost  of  this  reach  of  the 
coast  except  for  the  northern  and  southern  reaches  where  low 
waves  predominate.  Storm  surge  generally  increases  southward 
from  a  value  near  the  study  area  mean  in  the  north,  while  tidal 
range  remains  high  throughout.  From  Cape  Lookout  to  Bogue  Inlet 
sediment  size  is  variable  but  averages  about  the  mean  for  the 
mid-Atlantic  study  region.  Between  Bogue  Inlet  and  central  Ashe 
Island  sediments  are  fine.  From  central  Ashe  Island  to  New 
Topsail  Inlet  the  sediments  are  coarse.  From  New  Topsail  Inlet 
to  Fort  Fisher  beach  sediments  are  fine.  From  Fort  Fisher  to 
Cape  Fear  sediments  average  about  the  mean.  From  Cape  Lookout  to 
Beaufort  Inlet  the  bar  pattern  is  one  bar  or  no  bars.  The 
remainder  of  the  reach  is  dominated  by  a  one-bar  system  with 
occasional  small  areas  of  alternating  one-  and  two-bar  systems. 
Overwash  penetration  distance  is  low  over  most  of  the  region 
except  for  the  northern  and  southern  ends.  Host  of  this  section 
of  the  coast  is  relatively  stable  or  very  slightly  eroding. 

Between  Cape  Fear  and  the  South  Carolina-North  Carolina 
border  the  coast  strikes  easterly  and  has  a  moderate  offshore 
slope.  Island  width  and  lagoon  width  are  narrow  along  this  reach 
that  contains  many  inlets.  The  frequency  of  large  waves  is 
moderate,  storm  surge  is  high,  and  tidal  range  is  the  highest  of 
any  reach  along  the  mid-Atlantic  coast.  Beach  sediment  size  is 
variable  but  averages  near  the  mean.  The  dominant  bar  pattern  is 


a  one-bar  system.  Overwash  penetration  is  average  and  roost  of 
the  coastline  is  stable  or  eroding  at  a  slow  rate. 

Geographic  Subregions 

When  Principal  Components  Analysis  is  run  on  smaller 
subregional  data  sets  the  degree  of  organization  in  the  data 
structure  seems  to  be  greater.  This  is  evidenced  by  the  fact 
that  the  first  four  eigenvectors  typically  explain  more  than  70% 
of  the  variance  in  the  subregional  data  sets  compared  to  about 
60%  for  the  data  set  using  the  entire  coastal  reach.  In 
addition,  the  first  eigenvector  typically  accounts  for  about  40% 
of  the  variance  for  subregional  data  sets  compared  to  values  of 
about  20%  for  the  entire  coastal  data  set.  In  spite  of  the 
greater  degree  of  organization  evident  at  the  local  scale  there 
are  few  significant  differences  with  the  models  predicted  by  the 
analyses  of  the  entire  coastal  data  set  for  these  subregions. 
Minor  differences  occur  due  to  the  greater  resolution 
capabilities  of  the  subregional  analyses  while  some  of  these 
higher  resolution  changes  may  become  masked  by  the  entire  coastal 
data  set  at  the  larger  regional  scale.  We  will  address  only 
those  cases  where  significant  differences  are  visible  in  the 
analyses  of  the  subregional  data  at  the  local  scale,  and  the 
discussion  will  proceed  from  north  to  south. 

No  significant  differences  occurred  between  Cape  Henlopen 
and  Chincoteague  Inlet  (subregion  A).  This  suggests  that  the 
spatial  variance  in  the  coastal  geomorphic  and  process  parameters 
studied  are  organized  on  a  sufficiently  large  scale  to  be 


extracted  from  analysis  of  the  entire  coastal  data  set.  Between 
Chincoteague  Inlet  and  middle  Parramore  Island  (northern  half  of 
subregion  B) ,  the  islands  have  a  somewhat  higher  profile  and 
lagoons  are  slightly  wider  than  predicted  from  analyses  of  the 
entire  coastal  data  set.  The  remainder  of  subregion  B  from 
Parramore  Island  to  the  Chesapeake  is  not  significantly  different 
from  the  model  developed  from  the  entire  data  set.  In  subregion 
C ,  between  the  Chesapeake  Bay  and  Cape  Hatteras,  there  are 
several  areas  where  adjustments  to  the  large-scale  model  are 
suggested  from  analyses  of  the  subregions.  In  the  reach  between 
the  Chesapeake  Bay  and  False  Cape.  Virginia,  the  subregion  data 
suggest  there  are  fewer  offshore  bars  and  narrower  islands. 
Between  False  Cape  and  the  Virginia/North  Carolina  border  the 
subregional  data  suggest  higher  profile  islands,  steeper  offshore 
slopes,  and  slightly  coarser  sediments  than  indicated  by  the 
entire  coastal  data  set. 

Between  Oregon  Inlet  and  Rodanthe  the  subregional  analysis 
suggests  that  a  higher  profile  island  is  more  characteristic, 
high  waves  are  more  frequent,  slopes  are  slightly  steeper,  and 
the  lagoons  are  even  wider  than  indicated  by  the  entire  coastal 
data  set.  Between  Rodanthe  and  Cape  Hatteras  the  subregional 
analysis  suggests  that  beach  sediments  are  slightly  finer, 
lagoons  are  wider,  and  fewer  bars  occur  offshore  than  indicated 
by  the  entire  coastal  data  set. 

In  subregion  D,  between  Cape  Hatteras  and  central  Ocracoke 
Island,  the  islands  have  higher  profiles,  coarser  sediments,  and 
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lower  overwash  penetration  distance  than  indicated  by  the 
analyses  of  the  entire  coastal  data  set.  In  subregion  E»  the 
subregional  analysis  for  the  area  between  Cape  Lookout  and 
Beaufort  Inlet  suggests  that  higher  waves  occur  and  that  the 
southern  half  of  the  Cape  Lookout  to  Beaufort  Inlet  region  has 
higher  profile  islands  while  the  northern  half  has  lower  profile 
islands  than  indicated  by  the  analysis  of  the  entire  coastal  data 
set.  Between  Bogue  Inlet  and  central  Ashe  Island  the  subregion 
data  suggest  that  offshore  slope  is  slightly  steeper  and  that 
beach  sediments  are  slightly  coarser  than  indicated  by  the  entire 
coastal  data  set.  Between  Fort  Fisher  and  New  Topsail  Inlet  the 
subregion  model  predicts  that  there  are  coarser  sediments  and  a 
lower  frequency  of  large  waves  than  indicated  by  the  model  bas-ed 
on  the  entire  coastal  data  set. 


SUMMARY  AMD  CLASSIFICATION 


Coastal  classification  models  have  been  proposed  by 
numerous  investigators  over  the  past  century.  Most  of  these 
investigators  attempted  to  develop  a  unified  classification 
model  to  encompass  coasts  worldwide  (Davis  1912,  Johnson  1919, 
Cotton  1952,  Valentin  1952,  McGill  1958,  Shepard  1963,  Davies 
1964,  Inman  and  Nordstrom  1971.  Dolan  et  al.  1972).  Most  of 
these  models  were  based  on  the  vertical  or  horizontal  flux  of 
the  shoreline.  Davis  (1912)  proposed  a  variation  on  his  'cycle 
of  erosion'  to  include  coastlines  and  their  relation  to  uplift 
and  erosion.  The  model  proposed  by  Inman  and  Nordstrom  (1971) 
is  solely  based  on  the  tectonic  environment  of  world  coasts. 
Cotton  (1952)  made  the  first  level  of  classification  on  tectonic 
bases  and  then  subdivided  these  according  to  active  geomorphic 
processes.  Other  classifications  focused  upon  the  relative 
balance  between  transgressive  and  regressive  shorelines  (Johnson 
1919,  Valentin  1952)  . 

Some  models  placed  the  major  emphasis  on  geomorphic 
processes  important  in  the  genesis  of  coastal  sediments  (Shepard 
1963.  Inman  and  Nordstrom  1971,  Dolan  et  al.  1972).  Davies 
(1964)  classified  coasts  according  to  the  marine  process 
environment,  that  is,  by  wave  and  tidal  regime.  Kearns  (1974) 
classified  sandy,  coastal-plain  coasts  according  to:  the 
presence  of  barrier  islands,  the  vertical  activity  of  the 
coastal  region,  horizontal  shoreline  dynamics,  and  relative 


sediment  balance.  These  schemes  have  little  regional  or  local 
applicability  (i.e.,  for  the  mid-Atlantic  coast)  because  they 
were  designed  to  accommodate  the  varied  tectonic  and  geomorphic 
environments  of  world  coastlines#  and#  hence#  they  are  overly 
general  for  high-resolution  classification. 

Several  of  the  world  classification  schemes  were  based  on 
regional  observations  and  have  a  marked  regional  overprint. 
Tanner  (1960)  used  observations  of  the  Florida  coast  to  develop 
a  coastal  classification  based  on  the  lateral  stability  of  the 
shoreline.  Price's  (1959)  model  was  based  on  observations  made 
on  Gulf  of  Mexico  shorelines.  It  focuses  on  the  role  of  marine 
processes  in  shaping  coastlines  and  redistributing  sediments. 

No  detailed  classification  models  have  been  proposed  for 
the  mid-Atlantic  barrier  coast.  In  fact,  few  classification 
models  focus  on  barrier  coastlines.  The  regionalized 
classification  models  of  Price  (1959)  and  Tanner  (I960),  the 
hierarchical  model  of  Dolan  et  al.  (1972)#  and  the  sandy# 
coastal  plain  model  of  Kearns  (1974)  relate  better  to  the 
mid-Atlantic  coastline  than  the  general  world  models#  but  still 
can  not  provide  an  effective  system  for  classifying  the 
geomorphic  variations  occurring  along  that  coast. 

The  mid-Atlantic  study  reach  is  a  microtidal#  transgressive 
barrier  coast  with  limited  to  moderate  sediment  supply. 
Analysis  of  15  process  and  geomorphic  variables  (Table  4) 
provides  the  quantitative  basis  for  a  process-response  model. 
Our  plan  was  to  divide  the  coast  into  broad  regional 
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compartments  of  similar  geomorphic  and  process  attributes. 
Figure  7  summarizes  the  results  of  the  quantitative  analyses 
(principal  component  analysis)  with  the  24  distinct  coastal 
types  recognized  along  the  mid-Atlantic  coast.  It  is  important 
to  note  that  the  magnitudes  of  the  variables  are  represented 
relative  to  the  range  of  values  occurring  along  the  mid-Atlantic 
coast,  and,  therefore  are  not  directly  applicable  to  areas 
outside  the  800-km  study  area.  Many  of  the  distinctions  between 
high  and  low  magnitude  of  the  variables  are  necessarily  small 
because  of  the  geomorphic  similarity  within  this  microtidal 
study  reach.  Individual  attributes  display  apparent  spatial 
organization  at  various  scales  (See  Figs.  3,4,5  and  Tables  5,6), 
but  little  organization  appears  in  the  interrelationships 
between  all  of  the  parameters.  Clearly,  a  unified 
classification  must  utilize  fewer  attributes  if  a  smaller  number 
of  final  classification  categories  is  desired.  The 
relationships  between  two  or  three  variables  exhibit  greater 
organization  (for  example,  see  Resio  et  al.  1977,  Vincent  et  al. 
1976,  Dolan  et  al.  1977.  and  Dolan  et  al.  1979). 

By  a  first  approximation,  the  mid-Atlantic  coast  can  be 
divided  into  major  geomorphic  types  based  on  large-scale 
morphology:  1)  mainland  coasts  and  attached  barriers;  2)  long, 
continuous  barriers  (greater  than  25  km  in  length)  with  few 
inlets;  and  3)  short,  discontinuous  barriers  with  frequent 
inlets.  Figure  7  illustrates  the  systematic  recurrence  of  this 
morphologic  sequence  along  the  coast.  Along  the  open  coast 
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immediately  south  of  the  Delaware  Bay  and  the  Chesapeake  Bay, 
the  barrier  beaches  are  welded  to  the  mainland  (Fig.  7,  segments 
18  and  24).  South  of  these  areas,  long,  continuous  barriers 
dominate  the  coastline  (Fig.  7,  segments  6-17  and  21-23). 
Finally,  short,  segmented  barriers  occur  through  the  remainder 
of  each  reach  (Fig.  7,  segments  1-5  and  19-20).  By  deleting 
some  of  the  variables  from  consideration  (i.e.,  tidal  range, 
storm  surge,  and  overwash  penetration  distance)  the  coast  can  be 
subdivided  into  eight  regions  of  similar  geomorphic  attributes. 
The  first  of  these,  the  northernmost  reach  of  the  study  area 
between  Cape  Henlopen  and  southern  Assateague  Island  (Fig.  7, 
segments  22-24),  is  composed  of  mainland  coast,  attached 
barriers,  and  long,  barriers.  This  first  segment  is 
characterized  by  steep  offshore  slopes,  coarse-grained 
sediments,  zero  or  one  offshore  bar,  high  island  profiles,  high 
wave  frequency,  a  slowly  eroding  coastline,  and  moderately  wide 
lagoons  and  islands.  The  second  geomorphic  segment  is  between 
southern  Assateague  Island  and  the  Chesapeake  Bay  (Fig.  7, 
segments  19-21).  This  region  is  characterized  by  short, 
discontinuous  barriers,  very  gentle  offshore  slopes, 
fine-grained  sediments,  one  bar,  low  island  profiles,  low  wave 
frequency,  rapidly-eroding  coastlines,  wide  islands,  and 
moderately-wide  lagoons.  The  third  geomorphic  segment  occupies 
the  reach  between  the  Chesapeake  Bay  and  Nags  Head  (Fig.  7, 
segments  14-18) .  This  segment  begins  at  the  north  as  a 
mainland-attached  barrier  beach  and  becomes  a  long,  continuous 


barrier  to  the  south.  Geomorphic  attributes  of  this  segment 
include  gentle  to  moderate  offshore  slopes,  coarse-  to 
medium-grained  sediments,  temporally  variable  one-  or  two-bar 
systems,  high  island  profiles,  low  wave  frequency, 
slowly-eroding  coastlines,  moderately-wide  lagoons  and  wide 
islands.  The  fourth  geomorphic  segment  of  the  study  area  lies 
between  Nags  Head  and  Rodanthe  (Fig.  7,  segments  11-13).  This 
segment  is  characterized  by  long  barriers  with  steep  offshore 
slopes  (except  for  the  area  near  Oregon  Inlet  where  slopes  are 
gentle),  moderately  coarse-grained  sediments,  variable  one-  or 
two-bar  systems,  low  island  profiles,  moderate  wave  frequency, 
rapidly-eroding  coastlines,  and  wide  islands  and  lagoons.  The 
fifth  geomorphic  segment  of  the  mid-Atlantic  coast  is  a  reach  of 
long  barriers  between  Rodanthe  and  Cape  Lookout  (Fig.  7, 
segments  8-10)  .  This  reach  is  characterized  by  moderate  to 
steep  offshore  slopes,  medium-  to  coarse-grained  sediments, 
temporally  variable  one-  or  two-bar  systems,  low  island 
profiles,  high  wave  frequency,  stable  or  slowly-eroding 
coastlines,  moderate  to  wide  islands,  and  wide  lagoons.  The 
sixth  geomorphic  reach  of  the  mid-Atlantic  coast  is  between  Cape 
Lookout  and  Beaufort  Inlet  and  is  the  small  area  sheltered  from 
wave  activity  by  Cape  Lookout  (Fig.  7,  segment  7).  This  reach 
is  characterized  by  moder ately-gentle  offshore  slopes, 
medium-grained  sediments,  no  bars,  low  island  profiles,  low  wave 
frequency,  slowly-eroding  coastlines,  and  moderately-wide 
lagoons  and  islands.  The  seventh  geomorphic  segment  of  the 
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profile  islands.  Following  stabilization,  island-wide  overwash 
ceased,  islands  became  high  profile  barriers,  and  beaches 
narrowed  and  steepened  (Dolan  1972).  Coastal  engineering  and 
shore  protection  projects  have  altered  sediment  transport  along 
the  barrier  beaches  in  other  ways.  In  particular,  sand 
entrapment  structures  such  as  jetties  and  groins  can  cause 
serious  depletion  of  longshore  sand  supply  to  downdrift  beaches 
(for  example,  see  Ocean  City  Inlet  jetties).  Beach  nourishment 
projects  provide  excess  sand  to  selected  beaches.  In  these 
ways,  man's  activities  on  barrier  islands  can  significantly 
alter  natural  rates  of  shoreline  erosion  and  accretion. 
Therefore,  it  is  becoming  increasingly  difficult  along  populated 
beaches,  such  as  the  mid-Atlantic  barriers,  to  separate  natural 
response  from  man-induced  response.  This  rapidly  adds 
complexity  to  classification  schemes. 

Our  data  from  the  microtidal  mid-Atlantic  beaches  should 
provide  a  solid  data  base  for  local  and  regional  surveys  and 
coastal  zone  planning  along  the  study  reach.  In  addition,  it 
should  serve  as  a  good  basis  for  future  classifications  schemes 
in  other  environments  of  varied  tidal  energy,  wave  climate,  and 
tectonic  settings;  as  well  aB  for  comparative  studies  of  other 
microtidal  coastlines. 


mid-Atlantic  coast  is  the  reach  between  Beaufort  Inlet  and  Mason 
Inlet  (Fig.  7,  segments  3-6).  The  northern  part  of  this  reach 
is  characterized  by  long  barriers,  while  the  southern  part  is 
composed  of  short  barriers  with  frequent  inlets.  The  seventh 
segment  is  characterized  by  moderately-steep  offshore  slopes 
(except  segment  5  on  Fig.  7),  variable  sediment  size,  one 
offshore  bar,  high  island  profiles,  high  wave  frequency,  stable 
coastlines,  and  narrow  lagoons  and  islands.  The  eighth,  and 
southernmost,  georoorphic  segment  in  the  study  area  occurs 
between  Mason  Inlet  and  the  North  Carolina  -  South  Carolina 
border.  This  segment  is  characterized  by  short  barriers, 
moderate  offshore  slopes,  one  offshore  bar,  medium-grained 
sediments,  low  island  profiles,  low  wave  frequency,  stable 
coastlines,  and  narrow  lagoons  and  islands. 

We  have  been  able  to  categorize  the  coast  into  a  small 
number  of  similar  geomorphic  reaches  by  using  subsets  of  the 
data  base,  however,  when  all  fifteen  parameters  are  used,  the 
number  of  geomorphic  classes  increases  to  twenty-four.  Our 
analyses  have  delineated  several  reasons  why  our  attempts  to 
develop  a  simple  classification  with  fewer  than  eight  classes 
have  failed.  First,  there  does  not  appear  to  be  a  clear 
coupling  of  process  and  response  variables.  In  many  geomorphic 
systems  this  problem  of  equifinality  in  determining  the 
geomorphic  product  is  quite  common,  that  is,  there  are  many 
process  routes  to  the  same  morphometric  form  (response).  Our 
data  are  capable  of  delineating  response  on  a  1-km  resolution 


but  can  not  designate  direct  associations  with  process 
(particularly  on  the  same  scales). 

A  second  problem  with  classifying  the  coast  is  the 
disparity  in  resolution  scale  between  process  and  geomorphic 
response  variables.  In  order  to  explain  the  high-resolution 
variations  in  geomorphic  variables  measured  from  maps  and  aerial 
photographs,  process  variables  must  be  measured  on  a  much  higher 
resolution  scale  than  was  done  in  this  study.  Such  data  would 
require  a  large  investment  in  resources,  well  beyond  the  scope 
of  our  study,  where  process  and  response  variables  could  be 
measured  at  similar  scales. 

Third,  a  major  cause  of  difficulty  in  establishing  a 
process-response  classification  is  the  inability  to  distinguish 
between  morphometric  features  caused  by  modern  processes  vs. 
those  associated  with  relict  processes.  For  example,  offshore 
slope,  island  width,  lagoon  width,  and  sediment  size,  are 
probably  associated  with  relict  phenomena.  On  the  other  hand, 
overwash  penetration  distance,  the  rate  of  shoreline  change,  and 
bar  number  probably  are  in  local  equilibrium  with  modern 
processes. 

Fourth,  man's  interference  with,  and  manipulation  of, 
natural  barrier  island  processes  also  has  had  considerable 
impact  on  island  morphology.  Among  the  most  notable  of  man's 
efforts  to  change  barrier  morphology  are  the  massive  dune 
stabilization  projects  of  the  past  50  years.  Prior  to  dune 
stabilization,  many  areas  of  the  mid-Atlantic  barriers  were  low 
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profile  islands.  Following  stabilization,  island-wide  overwash 
ceased,  islands  became  high  profile  barriers,  and  beaches 
narrowed  and  steepened  (Dolan  1972).  Coastal  engineering  and 
shore  protection  projects  have  altered  sediment  transport  along 
the  barrier  beaches  in  other  ways.  In  particular,  sand 
entrapment  structures  such  as  jetties  and  groins  can  cause 
serious  depletion  of  longshore  sand  supply  to  downdrift  beaches 
(for  example,  see  Ocean  City  Inlet  jetties).  Beach  nourishment 
projects  provide  excess  sand  to  selected  beaches.  In  these 
ways,  man's  activities  on  barrier  islands  can  significantly 
alter  natural  rates  of  shoreline  erosion  and  accretion. 
Therefore,  it  is  becoming  increasingly  difficult  along  populated 
beaches,  such  as  the  mid-Atlantic  barriers,  to  separate  natural 
response  from  man-induced  response.  This  rapidly  adds 
complexity  to  classification  schemes. 

Our  data  from  the  microtidal  mid-Atlantic  beaches  should 
provide  a  solid  data  base  for  local  and  regional  surveys  and 
coastal  zone  planning  along  the  study  reach.  In  addition,  it 
should  serve  as  a  good  basis  for  future  classifications  schemes 
in  other  environments  of  varied  tidal  energy,  wave  climate,  and 
tectonic  settings;  as  well  as  for  comparative  studies  of  other 
microtidal  coastlines. 
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